NASA/CR— 2003-212351 



Quantification of Energy Release in 
Composite Structures 

Levon Minnetyan 

Clarkson University, Potsdam, New York 


June 2003 



The NASA STI Program Office ... in Profile 


Since its founding, NASA has been dedicated to 
the advancement of aeronautics and space 
science. The NASA Scientific and Technical 
Information (STI) Program Office plays a key part 
in helping NASA maintain this important role. 

The NASA STI Program Office is operated by 
Langley Research Center, the Lead Center for 
NASA's scientific and technical information. The 
NASA STI Program Office provides access to the 
NASA STI Database, the largest collection of 
aeronautical and space science STI in the world. 
The Program Office is also NASA's institutional 
mechanism for disseminating the results of its 
research and development activities. These results 
are published by NASA in the NASA STI Report 
Series, which includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant 
phase of research that present the results of 
NASA programs and include extensive data 
or theoretical analysis. Includes compilations 
of significant scientific and technical data and 
information deemed to be of continuing 
reference value. NASA's counterpart of peer- 
reviewed formal professional papers but 
has less stringent limitations on manuscript 
length and extent of graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies 
that contain minimal annotation. Does not 
contain extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by 
NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, 
often concerned with subjects having 
substantial public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific 
and technical material pertinent to NASA's 
mission. 

Specialized services that complement the STI 
Program Office's diverse offerings include 
creating custom thesauri, building customized 
databases, organizing and publishing research 
results . . . even providing videos. 

For more information about the NASA STI 
Program Office, see the following: 

• Access the NASA STI Program Home Page 
at http:/ hvww.sti.nasa.gov 

• E-mail your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA Access 
Help Desk at 301-621-0134 

• Telephone the NASA Access Help Desk at 
301-621-0390 

• Write to: 

NASA Access Help Desk 

NASA Center for AeroSpace Information 

7121 Standard Drive 

Hanover, MD 21076 



NASA/CR— 2003-212351 



Quantification of Energy Release in 
Composite Structures 

Levon Minnetyan 

Clarkson University, Potsdam, New York 


Prepared under Grant NAG3-2393 


National Aeronautics and 
Space Administration 


Glenn Research Center 


June 2003 



Acknowledgments 


This report is the documentation of work done at Clarkson University under NASA Grant NAG3-2393 (Clarkson 
University Research Project 375-32636 with Old Number 375-570 prior to 7/1/2002). Project dates were July 1, 2000 
to September 30, 2002 and total funding was $69,529. The project monitor was Dr. Christos C. Chamis, 

NASA Glenn Research Center. Many helpful discussions with Dr. Chamis on the implementation of 
this project are hereby acknowledged. Chapters I, II, and III are based on the graduate research of 
Mr. Hui Zhang who is a PhD student and who also helped with the assembly of this report. 

Chapters IV and V are based on the work of Ms. Prachi Suri who defended her 
MS degree at Clarkson in December 2002. 


This report is a formal draft or working 
paper, intended to solicit comments and 
ideas from a technical peer group. 


This report contains preliminary 
findings, subject to revision as 
analysis proceeds. 


The Propulsion and Power Program at 
NASA Glenn Research Center sponsored this work. 


Available from 


NASA Center for Aerospace Information 
7121 Standard Drive 
Hanover, MD 21076 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22100 


Available electronically at http: / /gltrs. grc.nasa.gov 


Table of Contents 


Page 

List of Tables vi 

List of Figures viii 

Chapter 1 . Energy Release of Adhesively Bonded Joints 1 

1.1 Introduction 1 

1 .2 Methodology 2 

1 .3 Method of Simulation 3 

1 .4 Simulation Results and Discussion 5 

1.5 Conclusions 8 

1.6 References 9 

Chapter 2. Energy Release Rates and Fracture Planes 17 

2.1 Introduction 17 

2.2 Methodology 20 

2.3 Microstress Level Damage Tracking 21 

2.4 Simulation of Braided Composite Notched Beam Specimens 22 

2.5 Conclusions 25 

2.6 References 26 

Chapter 3. Computational Prediction of Effective Elastic Constants in a Cross-Ply Laminate 

Under Uniaxial Loading 32 

3.1 Effect of 0° Constraint Plies on Stiffness Degradation in Cross-Ply Laminates 33 

3.2 Theoretical Model 35 

3.2.1 Lamina 1(90°) 35 

3.2.2 Lamina 2(0°) 38 

3.2.3 Strain Energy 41 

3.3 Effective Moduli 50 

3.4 Conclusions 51 

3.5 References 51 

Chapter 4. Calibration and Validation of Software Program ICAN/JAVA 53 

4.1 Objective 53 

4.2 Background 53 

4.3 Features and Enhancements of ICAN/JAVA 54 

4.4 Theory 55 

4.5 Model Description 58 

4.5.1 Durability Analysis 58 

4.5.2 Mfim Analysis 59 

4.6 Results and Discussion 60 

4.6.1 Cyclic in Plane Membrane Load 61 

4.6.2 Cyclic in Plane Bending Load 65 

4.6.3 Mfim Results 68 

4.7 Conclusion 70 

4.8 Future W ork 71 

4.9 References 71 


iii 


NASA/CR— 2003-212351 



Chapter 5. Prediction of First Ply Failure and Fracture in Composite Materials of Different Size 


and Geometry 74 

5.1 Objective 74 

5.2 Background 74 

5.3 Theory 75 

5.4 Model Description 77 

5.5 Model Assumptions 81 

5.6 Material Calibration 84 

5.7 Results and Discussion 88 

5.8 Conclusion 93 

5 . 9 Future W ork 93 

5.10 Material Properties 94 

5.1 1 References 98 


NASA/CR— 2003-212351 


IV 



List of Tables 


Page 

Table 1.1: AS-4 Fiber Properties 10 

Table 1.2: Epoxy Matrix Properties 10 

Table 2.1: X5260 Fiber Properties 27 

Table 2.2: G40-800 Matrix Properties 27 

Table 4. 1 : Durability Analysis Loads for the Test Cases 59 

Table 4.2: Loads Applied For MFIM Test 60 

Table 4.3: MFIM Input Options as Shown in IC AN/JAVA 60 

Table 4.4: Factor of Safety in the Longitudanal Direction For Cyclic in Plane 

Membrane Load 62 

Table 4.5: Factor of Safety in the Transverse Direction For Cyclic in Plane 

Membrane Load 63 

Table 4.6: Factor of Safety in the Shear Direction For Cyclic in Plane 

Membrane Load 64 

Table 4.7: Factor of Safety in the Longitudanal Direction for in Plane Bending Load 65 

Table 4.8: Factor of Safety in the Transverse Direction for in Plane Bending Load 66 

Table 4.9: Factor of Safety in the Shear Direction for in Plane Bending Load 67 

Table 5.1: Experimental and Simulated Data for Material Calibration 85 

Table 5.2: SGE1 -Fiber Properties 87 

Table 5.3: MHD1 -Matrix Properties 87 

Table 5.4: Comparison of Fiber Properties with Respect to Reference Material 88 

Table 5.5: Comparison of Matrix Properties with Respect to Reference Material 88 

Table 5.6: Coupon Simulation Data 88 

Table 5.7: Cylinder Simulation Data: SGES/MHDY Laminate 89 

Table 5.8: Cylinder Simulation Data: SGE1/MHD1 Laminate 90 


NASA/CR— 2003-212351 


v 



List of Figures 


Page 


Figure 1.1: Definitions of Regions for Ply Microstress Calculations 11 

Figure 1.2: Multi-Slice Unit Cell Subdivisions 12 

Figure 1 .2a: Unit Cell Square Array Concepts of Composite Micromechanics 11 

Figure 1.2b: Horizontal Slicing 12 

Figure 1.2c: Vertical Slicing 12 

Figure 1.3: Testing Setup for Adhesive Bonded Specimen 12 

Figure 1.4: Finite Element Model of Composite Specimen 12 

Figure 1.5: Displacement with Loading for Coupon 13 

Figure 1 .6: Displacement with Microstress Level Damage Energy Components 14 

Figure 1.6a: Displacement with SM1A(+) for Epoxy Resin Bonded Joints 13 

Figure 1.6b: Displacement with SM1A(+) for Graphite/Epoxy Prepreg Bonded Joint 13 

Figure 1 .6c: Displacement with SM2A(+) for Epoxy Resin Bonded Joints 13 

Figure 1.6d: Displacement with SM2A(+) for Graphite/Epoxy Prepreg Bonded Joints 13 

Figure 1.6e: Displacement with SM2B(-) for Epoxy Resin Bonded Joints 14 

Figure 1.6f: Displacement with SM2B(-) for Graphite/Epoxy Prepreg Bonded Joints 14 

Figure 1.6g: Displacement with SM13A for Epoxy Resin Bonded Joints 14 

Figure 1 .6h: Displacement with SM13A for Graphite/Epoxy Prepreg Bonded Joints 14 

Figure 1.7a: The Energy Envelope of Acoustic Emission Signals for Epoxy Resin 

Bonded Specimen 15 

Figure 1.7b: The Energy Envelope of Acoustic Emission Signals for Prepreg 

Bonded Specimen 15 

Figure 1.8: Damage with Damage Energy for Coupon 16 

Figure 2.1a: Configuration for E-Type Specimen 28 

Figure 2.1b: Configuration for IE-Type Specimen 28 

Figure 2.2a: Finite Element Model (Brick Element) for E-Type Specimen 29 

Figure 2.2b: Finite Element Model (Shell Element) for IE-Type Specimen 29 

Figure 2.3: Displacement with Loading for IE-Type and E-Type Coupons 30 

Figure 2.4: Damage Energy with Damage Volume for IE-Type and E-Type Coupons 30 

Figure 2.5: Loading with Damage Volume for IE-Type and E-Type Coupons 31 

Figure 2.6: Fracture Toughness with Crack Area for IE-Type and E-Type Coupons 3 1 

Figure 3.1a: Schematic of A [0790°] s Laminate with Transverse Cracks Under 

Uniaxial Loading 33 

Figure 3.1b: One-Quarter of Unit Damaged Cell with Transverse Cracking 34 

Figure 4. 1 : Composite Ply Lay-Up for Durability Analysis 58 

Figure 4.2: Composite Ply Lay-Up for Mfim Analysis 59 

Figure 4.3: Graph of Ml 1 for Cyclic in Plane Membrane Load 62 

Figure 4.4: Graph of M22 for Cyclic in Plane Membrane Load 63 

Figure 4.5: Graph of M12 for Cyclic in Plane Membrane Load 64 

Figure 4.6: Graph of Ml 1 for Cyclic in Plane Bending Load 66 

Figure 4.7: Graph of M22 for Cyclic in Plane Bending Load 67 

Figure 4.8: Graph of M12 for Cyclic in Plane Bending Load 68 

Figure 4.9: Graph of Ml 1 for MFIM 69 

Figure 4.10: Graph ofM22 for MFIM 69 


NASA/CR— 2003-212351 


vi 



Figure 4.1 1: Graph of M 12 for MFIM 70 

Figure 5.1: Simulation Cycle of Progressive Fracture in Composite Laminates 

and Structures Via CODSTRAN 75 

Figure 5.2: CODSTRAN Load Increment 76 

Figure 5.3: Overall CODSTRAN Simulations 77 

Figure 5.4: E-glass/MY750 [±55] s Coupon 78 

Figure 5.5: E-glass/MY750 [±55] s Cylindrical Tube 79 

Figure 5.6: Twelve Ply Schedules to Simulate Reinforcement for Cylinder 80 

Figure 5.7: Variable Thickness of the Cylinder Due to Different Ply-Schedules 80 

Figure 5.8: Boundary Conditions for the Coupon 81 

Figure 5.9: E-glass/MY750 Tube Geometry 82 

Figure 5.10: Direction of Force and Pressure Applied to the Cylinder Specimen 83 

Figure 5.11: Biaxial Envelope of 0° GRP Lamina Under Combined Cx and Gy Stresses 

for E-glass/MY750 Epoxy 84 

Figure 5.12: Comparison of Biaxial Envelope of 0° GRP Lamina Under Combined Ox and Gy 

Stresses for Experimental and Simulated Material Properties 86 

Figure 5.13: Biaxial Failure Stresses for [±55] s E-glass/My750 Laminate 

Using CODSTRAN 91 

Figure 5.14: Biaxial Failure Stresses for [±55] s E-glass/My750 Laminate 

Obtained Experimentally 92 


vii 


NASA/CR— 2003-212351 



Biographies 


Mr. Hui Zhang is a PhD student in the Department of Civil and Environmental Engineering at 
Clarkson University. He received his M.S. degree in Structural Engineering from the Hefei 
University of Technology in 1998. His research area is in computational simulation of 
progressive fracture of composite structures. 

Dr. Levon Minnetyan is a professor of structural engineering at Clarkson University in the 
Department of Civil and Environmental Engineering. His primary research activities are directed 
toward the assessment of progressive damage and fracture in structures made of composite 
materials such as graphite/epoxy laminated, woven, and braided composites. 


viii 


NASA/CR— 2003-2 1235 1 



Chapter 1 

Energy Release of Adhesively Bonded Joints 


In this chapter, the adhesive bond strength of lap-jointed graphite/aluminum composites is 
examined by computational simulation. Computed micro-stress level energy release rates are 
used to identify the damage mechanisms associated with the corresponding acoustic emission 
(AE) signals. Computed damage regions are similarly correlated with ultrasonically scanned 
damage regions. Results show that computational simulation can be used with suitable NDE 
methods for credible in-service monitoring of composites. 

Keywords: Computational simulation, Energy release rate, Nondestructive Evaluation (NDE) 

1.1 Introduction 

For effective structural health monitoring, it is important to quantify damage tolerance of a 
candidate structure. Since continuous fiber composites are able to arrest cracks and prevent self- 
similar crack propagation, composite structures have received a great deal of consideration for 
design with emphasis on damage tolerance. However, a number of design parameters such as 
fiber orientation patterns, choices of constituent material combinations, ply drops and 
hybridization, result in complex design options for composite structures. Thus, it is necessary to 
evaluate damage initiation in a composite structure and its fracture propagation characteristics 
for achieving a rational damage tolerant design. 

Compared with homogeneous materials, damage initiation and progression characteristics of 
fiber composites are much more complicated. Composite structures often contain some pre- 
existing or induced flaws in matrix and fibers after fabrication of composites. At lower stresses, 
matrix is likely to be cracked because of flaw-induced stress concentrations and cause the matrix 
flaws to propagate across the composite. With the use of established material modeling and finite 
element models, and considering the influence of local defects, through-the-thickness cracks and 
residual stresses, computational simulations have made it possible to evaluate the details of 
progressive damage and fracture in composite structures. In a computational simulation, damage 
evolution quantifier such as the damage volume, exhausted damage energy, and the damage 
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energy release rate (DERR) are used to quantify the structural damage tolerance at different 
stages of degradation. Low DERR levels usually indicate that degradation takes place with minor 
resistance by the structure. Structural resistance to damage propagation is often dependent on 
structural geometry and boundary conditions as well as the applied loading and the state of 
stress. 

In certain cases such as the room temperature behavior of composites designed for high- 
temperature applications, internal damage initiated as microcracks in the matrix become enlarged 
to be externally visible. Thus, matrix cracking and its effect on damage propagation/damage 
tolerance need be evaluated. Some simulations [1,2] have been successful in predicting damage 
tolerance and failure load of composite structure by considering ply stresses and the 
corresponding stress limits for matrix crack growth. In this report, lap-jointed composite 
specimens subjected to uniaxial tension are investigated. Damage initiation, growth, 
accumulation, and propagation to fracture are studied. Since the complete evaluation of ply and 
subply level damage/fracture processes is the fundamental premise of computational simulation, 
a microstress level damage index is added for the identification and tracking of subply level 
damage processes. Computed damage regions are similarly correlated with ultrasonically 
scanned damage regions. Simulation is validated by comparison with test data from acoustic 
ultrasonic (AU) testing [7]. Results show that computational simulation can be used with suitable 
NDE methods for credible in-service monitoring of composites. 

1.2 Methodology 

Computational simulation is implemented via integrating three modules: (1) composite 
mechanics, (2) finite element analysis, and (3) damage progression tracking. The composite 
mechanics module (Murthy and Chamis 1986) is designed to analyze fiber composite structures 
with an updated composite mechanics theory. Its main function is to calculate ply and composite 
properties of laminates from the fiber and matrix constituent characteristics and the composite 
layup. Additionally, it determines the composite structural response and ply stresses from the 
FEM analysis results. In simulation, the composite mechanics module is called before and after 
each finite element analysis. 
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The finite element analysis module is able to process linear and nonlinear static and dynamic 
analysis. Four-node anisotropic thick shell elements are usually used to model laminated 
composites (Nakazawa et al 1987). The finite element analysis module accepts laminate 
properties from the composite mechanics module and performs the structural analysis at each 
load increment. After structural analysis, the computed generalized node stress resultants and 
deformations are provided to the composite mechanics module. The composite mechanics 
module computes the developed ply stresses for each ply and checks for ply failure modes at 
each node. Failure criteria applied to detect ply failures are based on the maximum stress and 
modified distortion energy (MDE) criteria for combined stress effects (Murthy and Chamis 
1986). 

The overall evaluation of composite structural durability is carried out in the damage 
progression module (Minnetyan et al 1990) that keeps track of composite degradation for the 
entire structure. The damage progression module relies on the composite mechanics module for 
composite micromechanics, macromechanics and laminate analysis, and calls the finite element 
analysis module for global structural analysis. If excessive damage is detected, the incremental 
loads are reduced and the analysis is restarted from the previous equilibrium stage. Otherwise, if 
the increment of loads is acceptable, another finite element analysis is perfonned but the 
constitutive properties and the finite element mesh are updated to account for the damage and 
deformations from the last simulation. Simulation is stopped when global structural fracture is 
predicted. 

1.3 Method of Simulation 

The matrix in orthotropic composite plies is divided into two parts: regions A and B. Region 
A represents the area in which stress concentrations induced by the interaction of matrix and 
fiber do not create any effect in matrix. Region B represents the interaction zone between fiber 
and matrix. Figure 1.1 shows the details of regional subdivision in transverse and normal 
directions of a composite ply with square packing. 

Considering the behavior of longitudinal stress <r m , transverse stress <J I22 , in-plane shear 

stress <7 n2 , out-of plane shear stress <J l23 , temperature gradient AT,, and moisture M,, Murthy 
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and Chamis (1986) present the complete set of equations for evaluating ply microstresses in 
regions A and B. For example, ply microstresses due to <7 in are given by: 


V,nU =(E m /E m )(T, n 
o /ii = \ E f\\ ! E m p m 

= O'ils = {v m ~ v m\ E ,n I E ,u)° in 


(7^ = <7^ = (7^ = (7^ = — 

^ mil ^ m33 u f 22 ^ f 33 


1 Jky 


1 11 


in which a mU is the matrix longitudinal stress, <T /n is the fiber longitudinal stress, is the 
matrix transverse stress in region A, and cr$ 2 is the matrix transverse stress in region B. If the 

ply is subject to combined stresses, its microstresses are obtained by simply superimposing 
results of all corresponding stress components. Ply transverse fractures usually begin in region B 
due to the elevated stress levels from stress concentration. Microstress level damage tracking is 
able to quantify the type of damage in the matrix by comparison of microstresses with 
constituent stress limits. A microstress damage index is defined as a binary number with 14 bits 
in the damage progression module. The 14-bit number corresponds to the following components: 
(1) SM1A(+) longitudinal stress in region A, positive; (2) SMIA(-) longitudinal stress in region 
A, negative; (3) SM2A(+) transverse stress in region A, positive; (4) SM2A(-) transverse stress 
in region A, negative; (5) SM2B(+) transverse stress in region B, positive; (6) SM2B(-) 
transverse stress in region B, negative; (7) SM3A(+) normal stress in region A, positive; (8) 
SM3A(-) normal stress in region A, negative; (9) SM3B(+) normal stress in region B, positive; 
(10) SM3B(-) nonnal stress in region B, negative; (11) SM12A in-plane shear stress in region A; 
(12) SM12B in-plane shear stress in region B; (13) SM13A out-of-plane shear stress in region A; 
(14) SM23A out-of-plane shear stress in region A. When the binary bit corresponding to 
SM2B(+) is set equal to 1, it indicates that region B is fractured with the transverse failure mode. 
In a subsequent stage transverse cracking will spread to region A. 

Mital (1993) presents a more refined slice-by-slice substructuring model for assessing ply 
failure modes at the microstress level. In his model, fiber is substructured into several slices and 
the micromechanics equations are applied at the slice level. Once the equivalent slice properties 
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are defined, the ply properties are obtained using classical lamination theory. By using maximum 
strength and combined stress criteria on the fiber, matrix, and interface microstresses, failure 
criteria are applied directly at the slice microstress level. With this model, the variation of 
interfacial bonding around the fiber circumference can be depicted more clearly. As a result, the 
local matrix cracking and fiber breaks can be monitored more closely in computational 
simulation. For ceramic matrix composites structural life prediction is quantifiable more 
accurately with the slice-by-slice ply substructuring model. The fiber substructuring and slice 
geometry are shown in Figure 1.2. 

1.4 Simulation Results and Discussion 

In this report, two unidirectional graphite/aluminum composite plates with single lap joint 
under uniaxial tension are used to demonstrate the use of microstress damage index in 
computational simulation. The specimen has a length of L = 76.2mm (3 in.), a width of W = 
19.05mm (0.75 in.) and a thickness of H = 1.27mm (0.05 in.). The fiber volume ratio is 
V f = 0.60 and the void volume ratio is V v = 0.01 . The specimens are bonded with adhesive to 

produce a bond area of 0.75"x0.75". The first adhesive type is an epoxy resin and the other is a 
graphite/epoxy prepreg tape. The specimen is assumed to be dry with zero moisture content. 
Additionally, two metal plates with dimensions of 2"x0.75" are also bonded to the specimens 
with epoxy resin. The fiber and matrix properties used for computational simulation are given in 
Table 1.1 and Table 1.2, and the configuration for adhesive bonded specimen is shown in Figure 
1.3. 

The finite element model as shown in Figure 1.4 has 451 nodes and 360 Mindlin type 
rectangular elements. To represent the test setup described by Quattlebaum (1997), nodes along 
the end of the specimen are restrained against translation to model the fixed boundary and nodes 
at the end of another specimen are constrained only moveable in the longitudinal direction. Axial 
tension load is applied at the moveable end. Moreover, nodes along the moveable end are tied by 
duplicate node constraints to enforce the uniform displacement of the clamped loaded edge. 
Computational simulation indicated a damage initiation load of 4555. 06N (1024 lbs) for the lap 
joint with epoxy resin. The damage initiation mode was by ply longitudinal tensile failures in the 
0° ply and the microstress damage indexes are activated for the SM1A(+), SM2A(+), SM2B(-), 
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SM3A(+), SM3B(-), SM12A, SM12B, SM13A, and SM23A microfailure modes. The locations 
of damage initiation were at the lap joint area. After damage initiation, longitudinal tensile 
failures spread to the interface between lap joint and composite plate as the applied load reached 
7584. 35N (1705 lbs). Then, longitudinal tensile failures continuously spread at the lap joint with 
the load increasing to 7691 N (1729 lbs). The microstress damage indexes were also activated for 
the SM1A(+), SM2A(+), SM2B(-), SM3A(+), SM3B(-), SM12A, SM12B, SMI 3 A, and SM23A 
microfailure modes. Thereafter, longitudinal tensile failures developed at the lap joint and the 
interface between lap joint and the composite plate as the load increased. The lap joint continued 
to fracture rapidly and broke under the 44.096kN (9913 lbs) loading. Computational simulation 
indicated that adhesive yielding resulted in failure of the lap joint. 

For the lap joint with graphite/epoxy prepreg tape as its adhesive, the damage initiation load 
was 36.45kN (8193 lbs) by ply longitudinal tensile failures in the 0° ply. The damage initiation 
took place at the edges of the lap joint adjacent to the composite plate. The microstress damage 
indices were activated for SM1A(+), SM2A(+), SM2B(-), SM3A(+), SM3B(-), SM12A, SM12B, 
SM13A, and SM23A microfailure modes. After damage initiation, longitudinal tensile failures 
spread to plies of the specimen as the applied load reached 42.52kN (9558 lbs). Subsequently, 
damage growth was limited mainly to the first two plies of interface between the edge of lap 
joint and the composite plate until the load was increased to 44.53kN (10010 lbs). At 44.75kN 
(10060 lbs), longitudinal tensile failures spread to plies around the edge of lap joint. The 
microstress damage indices were activated for SM1A(+), SM2A(+), SM2B(-), SM3A(+), 
SM3B(-), SM12A, SM12B, SM13A, and SM23A microfailure modes. With increasing load, 
fracture continued to develop at the edges of the lap joint and the specimen broke under the 
52.49kN (11800 lbs) loading. Similar to the simulation for the lap joint bonded with epoxy resin, 
adhesive yielding caused failure of bond joint. 

Figure 1.5 shows the relation between displacement and loading for lap joint with two 
different adhesives. It indicates the bonded strength of lap joint with graphite/epoxy prepreg tape 
is approximately 1.19 times higher than that of the epoxy bonded joint, a little lower than 
experimental results by Quattlebaum (1997). This is due to some difference in the material 
properties in simulation and experiment. On the other hand, the apparent linear relationship 
between the load and displacement depicted in Figure 1.5 is not able to reflect the presence of 
internal damage initiation and growth processes. 
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For epoxy and graphite/epoxy prepreg tape bonded specimens, Figure 1.6a and Figure 1.6b 
compare the change of microstress energy component SM1A(+) with increasing displacement. It 
shows that the peak of SM1A(+) for graphite/epoxy prepreg tape is quite later and higher than 
that of epoxy. This means cracking in region A for epoxy bonded joint is earlier than prepreg 
bonded joint. The difference in the damage energy amplitudes indicate that the prepreg bonded 
joint has a much higher energy of SM1A(+) damage. 

Figure 1.6c and Figure 1.6d show the microstress energy component SM2A(+) as a function 
of displacement. From the plots, it is obvious that prepreg specimen can withstand higher 
transverse tensile stress in large displacement than epoxy specimen. It reflects that there is 
greater resistance against crack initiation in matrix for prepreg specimen. 

Figure 1.6e and Figure 1.6f show the relation of microstress energy component SM2B(-) 
versus displacement. It is observed that the stress required to cause transverse debonding at the 
fiber-matrix interface for prepreg specimen is larger than that for epoxy specimen. The bond 
strength for epoxy specimen is not as strong as that for prepreg specimen. This is mainly due to 
the fiber bridging effects that impede crack propagation in the prepreg-bonded specimen. 

Figures 1.6g and Figure 1.6h plot microstress energy component SM13A versus 
displacement. In the graph, the shear stress (out of plane) in region A for prepreg bonded joint is 
higher than that for epoxy bonded joint. It shows that there exist bridging fibers in prepreg 
specimen. 

Recently Nondestructive Evaluation (NDE) methods have been received more attention. Raju 
(1996) studied the in-situ process of C-C composites with Acoustic-Ultrasonic method. 
Quattlebaum (1997) presented the acoustic activity data obtained from epoxy resin and prepreg 
specimens with the same configuration as the computational simulation. Direct comparison of 
the acoustic emission signals with computed microstress level damage energies was not possible 
because the computed damage energy values had much fewer points than the experimentally 
measured acoustic emission data. The apparent reason for the mismatch of refinement between 
experimental and computational values is that microstress level damage tracking is not sufficient 
to capture the acoustic emission data from NDE testing. During the proposed continuation of this 
project a nanoscale level damage energy tracking will be implemented using progressive 
refinement of the composite subregions via new software that have been implemented in the 
IC AN- Java code. IC AN- Java is a new generation of composite mechanics code that enables the 
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telescopic multi-scale subdivision of composites in both the thickness and the 3-D spatial sense. 
Integration of ICAN-Java refinement capability in the CODSTRAN progressive damage and 
fracture code will be necessary to render progressive decomposition of the computationally 
simulated damage energies to NDE signal precision. For correlation of current microstress level 
damage energy tracking with Quattlebaum’s test data, we used a curve fitting software 
(TableCurve2D v5 2000) to find the fitted curve that envelops acoustic signals. Figure 1.7a 
shows the envelope curve for epoxy resin bonded specimen and Figure 1.7b plots the envelope 
curve for prepreg bonded specimen. Experimental acoustic emission levels are marked with x as 
points through which curves are drawn in Figures 1.7a and 1.7b. The area under the envelope 
represents the total damage energy detected by the acoustic emission during the period 
monitored. Thus, we can correlate the microstress damage energy through computational 
simulation with the relative total damage energies represented by the envelope areas. The 
envelope area in Figure 1.7a is 1.22986e-3 and that in Figure 1.7b is 1.30145e-3. The ratio of the 
energies from Figure 1.7b to that of Figure 1.7a is 1.058. Comparatively, the ratio of simulation 
results for damage energy shown in Figure 1.8 is 1.168, which is 10 percent higher than the ratio 
computed from the NDE test results. 

1.5 Conclusions 

The failure patterns of the investigated fiber composite specimens and the available 
computational simulation results are found that: 

1. Microstress level damage tracking is able to evaluate the bond strength and monitor 
adhesive yielding. 

2. Microstress level computational simulation provides the details of damage initiation, 
growth, and subsequent fracture in composites. It represents a new approach for 
investigating damage mechanisms of composites. 

3. Computational simulation shows a good correlation with AU signals in the cumulative 
sense. 

4. More refined nanoscale damage energy tracking via progressive scale decomposition 
implemented in the ICAN-Java code will be necessary to enable direct comparison with 
NDE AU signals. 
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5. The demonstrated procedure is flexible and applicable to all types of constituent 
materials, structural geometry, and loading. Hybrid composites, as well as laminated, 
stitched, woven, and braided composites can be simulated. 

6. Computational simulation, with the use of established composite mechanics and finite 
element models, can be used to predict the influence of microstresses, as well as loading 
and material properties on the durability of composite structures. 


1.6 References 


1. H. Zhang, L. Minnetyan and C.C. Chamis (2001) “Microstress Level Damage Evolution in 
Composite Structures,” 46 th International SAMPE Symposium, Long Beach, CA 

2. L. Minnetyan, P.L.N. Murthy, and C.C. Chamis (1990) “Composite Structure Global 
Fracture Toughness via Computational Simulation,” Computers & Structures, Vol. 37, No. 
2, pp. 175- 180 

3. William H. Prosser (1996) “Advanced AE Techniques in Composite Materials Research”, 
Journal of Acoustic Emission . Vol. 14(3-4), pp.Sl-Sll 

4. P.L.N. Murthy and C.C. Chamis (1986) Integrated Composite Analyzer (ICAN): Users and 
Programmers Manual , NASA Technical Paper 2515 

5. S. Nakazawa, J.B. Dias, and M.S. Spiegel (1987) MHOST Users’ Manual , Prepared for 
NASA Lewis Research Center by MARC Analysis Research Corp., April 

6. S.K. Mital, P.L.N. Murthy, and C.C. Chamis (1993) “Ceramic Matrix Composites Properties/ 
Microstresses with Complete and Partial Interphase Bond,” NASA Technical Memorandum 
106136, 12 p. 

7. G.J. Quattlebaum and P.K. Raju, “Nondestructive Evaluation of Graphite/Aluminum 
Composites Using Acousto-Ultrasonics,” Proceedings of the ASME Noise Control and 
Acoustics Division, ASME 1997. 

8. P.K. Raju and U.K. Vaidya (1996), “Acoustic-Ultrasonic Evaluation of Carbon-Carbon 
Composites at Different Processing Stages,” Polymer Composites , Vol. 17, No.2, pp275-287. 

9. TableCurve2D v5 (2000), AISN Software. 


NASA/CR— 2003-212351 


9 



Table 1.1: AS-4 Fiber Properties 


Number of fibers per end = 10000 
Fiber diameter = 0.00508 mm (0.200E ' 3 in) 

Fiber Density = 4.04E' 7 Kg/m 3 (0.063 lb/in 3 ) 

Longitudinal normal modulus = 226.84 GPa (3.29E +7 psi) 

Transverse normal modulus = 13.72 GPa (0.199E +7 psi) 

Poisson's ratio (v 12 ) = 0.200 
Poisson's ratio ( v 23 ) = 0.250 

Shear modulus (G 12 ) = 13.79 GPa (0.20E +7 psi) 

Shear modulus (G 23 ) = 6.89 GPa (0.10E +7 psi) 

Longitudinal thermal expansion coefficient = 1.0E' 6 /°C (-0.55E' 6 /°F) 
Transverse thermal expansion coefficient = 1.0E' 6 /°C (0.56E' 6 /°F) 
Longitudinal heat conductivity = 301 kJ-m/hr/m 2 /°C (4.03 BTU-in/hr/in 2 /°F) 
Transverse heat conductivity = 30.1 kJ-m/hr/m 2 /°C (0.403 BTU-in/hr/in 2 /°F) 
Heat capacity = 0.712 kJ/kg/°C (0.17 BTU/lb/°F) 

Tensile strength = 3.72 GPa (540 ksi) 

Compressive strength = 3.35GPa (486 ksi) 


Table 1.2: Epoxy Matrix Properties 

Matrix density = 3.30E' 7 Kg/m 3 (0.0443 lb/in 3 ) 

Normal modulus = 3.45 GPa (500 ksi) 

Poisson's ratio = 0.35 

Coefficient of thermal expansion = 0.77E' 4 / °C (0.428E' 4 / °F) 

Heat conductivity = 0.648 kJ-m/hr/m 2 / °C (0.868E ' 2 BTU- in/hr/in 2 / °F) 
Heat capacity = 1.047 KJ/Kg/ °C (0.25 BTU/lb/ °F) 

Tensile strength = 68.99 MPa (10.0 ksi) 

Compressive strength = 241.59 MPa (35.0 ksi) 

Shear strength = 89.7 MPa (13.0 ksi) 

Allowable tensile strain = 0.02 
Allowable compressive strain = 0.05 
Allowable shear strain = 0.045 
Allowable torsional strain = 0.045 

Void conductivity = 16.8 J-m/hr/m 2 / °C (0.225 BTU-in/hr/in 2 / °F) 
Glass transition temperature = 216 °C (420°F) 
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Figure 1.2a: Unit Cell Square Array Concepts of Composite Micromechanics 
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Figure 1.2b: Horizontal Slicing Figure 1.2c: Vertical Slicing 


Figure 1.2: Multi-Slice Unit Cell Subdivisions 



Figure 1.3: Testing Setup for Adhesive Bonded Specimen 


Figure 1.4: Finite Element Model of Composite Specimen 
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Figure 1.5: Displacement with Loading for Coupon 
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Figure 1.6a: Displacement with SM1A(+) 
for Epoxy Resin Bonded Joints 



Figure 1.6b: Displacement with 
SM1A(+) for Graphite/Epoxy 
Prepreg Bonded Joints 
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Figure 1.6d: Displacement with 
SM2A(+) for Graphite/Epoxy 
Prepreg Bonded Joints 
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Figure 1.6f: Displacement with SM2B(-) 
for Graphite/Epoxy Prepreg Bonded Joints 



Figure 1.6g: Displacement with SM13A for 
Epoxy Resin Bonded Joints 
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Figure 1.6h: Displacement with SM13A for 
Graphite/Epoxy Prepreg Bonded Joints 


Figure 1.6: Displacement with Microstress Level Damage Energy Components 
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Figure 1.7a: The Energy Envelope of Acoustic Emission 
Signals for Epoxy Resin Bonded Specimen 
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Figure 1.7b: The Energy Envelope of Acoustic Emission 
Signals for Prepreg Bonded Specimen 
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Chapter 2 

Energy Release Rates and Fracture Planes 


Energy release rate G is usually suggested as a parameter for characterizing material 
toughness. Computational prediction of energy release rate is based on composite mechanics 
with micro-stress level damage assessment, finite element structural analysis and damage 
progression tracking modules. In this chapter, mode I interlaminar and intralaminar energy 
release rates G IC of composites are examined by computational simulation. Results show that 

computational simulation has a good predictive capability for monitoring damage progression in 
composites. Computational simulation enables assessment of the damage initiation and 
propagation loads. Computational simulation can be used prior to testing. Through simulation, 
sensitive parameters affecting critical values of the energy release rates are identified, which 
significantly enhance the accuracy and productivity of experiments. Simulation results are 
compared with test data. 

Keywords: composites, composite materials, energy release rates, notched beam specimen, 
computational simulation, progressive damage 

2.1 Introduction 

Critical components of a structure are required to remain safe and be able to function under 
loading after experiencing some damage. The cause of damage may be an accident, defect, or 
unexpected overloading. Damage tolerance of a structure is quantified by the residual strength, 
that is the additional load carrying ability after damage. Design considerations with regard to the 
durability of metallic and fiber composite structures require an a priori evaluation of damage 
initiation and fracture propagation mechanisms under expected loading and service 
environments. Concerns for safety and survivability of critical components require a 
quantification of the structural fracture resistance under loading. For a rational design process it 
is necessary to quantify the structural damage tolerance for a candidate design. The ability of 
designing composites with numerous possible fiber orientation patterns, choices of constituent 
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material combinations, ply drops and hybridizations, render a large number of possible design 
parameters that may be varied for an optimal design. The structural fracture process of a fiber 
composite depends on many parameters such as laminate configuration, fiber volume ratio, 
constituent stiffness/strength/hygro thermal parameters, stiffening system, and the fabrication 
process. Recent developments in computational simulation technology implemented in the GRC 
CODSTRAN (Composite Durability STRuctural ANalysis) code have made it possible to 
evaluate the details of progressive damage and fracture in metallic and composite structures. 
Computational simulation enables assessment of the damage initiation and propagation loads. 
The influence of local defects or flaws and effects of the fabrication process in terms of residual 
stresses are also taken into account. 

Computational simulation of 2-D and 3-D fiber composite structures have been successful in 
quantifying the load-displacement relationships, fracture paths, and ultimate strength of 
composite structures. Computational simulation perfonns a complete evaluation of laminated 
composite fracture via assessment of ply and subply level damage/fracture processes. The 
evaluation of composite test response can be made much more productive and informative via 
computational evaluation of energy release rates. Computational simulation can be used prior to 
testing of a metallic or fiber composite structure for evaluation of changes in the local and global 
critical stress intensity factors and strain energy release rates. Progressive damage mechanisms, 
damage locations/modes, and sensitive parameters affecting failure can be identified prior to 
testing, significantly enhancing the accuracy and productivity of an experimental program. 

An important feature of computational simulation is the assessment of damage stability or 
damage tolerance of a structure under loading. At any stage of damage progression, if there is a 
high level of structural resistance to damage progression under the service loading, the structure 
is stable with regard to fracture. The corresponding state of structural damage is referred to as 
stable damage. On the other hand, if damage progression does not encounter significant 
structural resistance, it corresponds to an unstable damage state. Unstable damage progression is 
characterized by very large increases in the amount of damage due to small increases in loading. 
Whereas during stable damage progression the amount of increase in damage is consistent with 
the increase in loading. 

Internal damage in composites is often initiated as cracking due to normal stresses transverse 
to fiber orientation. At the presence of stress concentrations or defects, initial damage may also 
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include fiber fracture. Further degradation is in the form of additional fiber fractures that usually 
lead to structural fracture. Because of the numerous possibilities with material combinations, 
composite geometry, fiber orientations, and loading conditions, it is essential to have an effective 
computational capability to predict the behavior of composite structures for any loading, 
geometry, composite material combinations, and boundary conditions. The predictions of 
damage initiation, growth, accumulation, and propagation to fracture are important in evaluating 
the load carrying capacity and reliability of composite structures. Quantification of the structural 
fracture resistance is also required to evaluate the durability/life of metallic and composite 
structures. 

Compared with homogeneous materials, damage initiation and progression characteristics of 
fiber composites are much more complicated. One of the most important problems is interface 
cracking that is known as delamination. Since composite structures often contain some pre- 
existing or induced flaws in matrix and fibers after fabrication of composites, delaminations 
frequently occur at lower stresses. For preventing the degradation in the strength and damage 
tolerance of a structure due to delaminations, it is necessary to quantify the effect of 
delaminations. 

In general, there are two fracture modes induced by delaminations; interlaminar and 
intralaminar fracture, that are observed in composites. Experimental research [1,2] has been 
extensively conducted on the interlaminar fracture mode. Correspondingly, a few works [3,4] 
attempted to investigate the behavior of intralaminar fracture mode. However, a number of 
design parameters such as fiber orientation patterns, choices of constituent material 
combinations, ply drops and hybridization, result in complex design options for composite 
structures in which intralaminar transverse tensile fractures usually precede more complex 
fracture modes. Thus, it is difficult for experiments to assess damage propagation in different 
situations. 

Recent developments in computational simulation technology [5,6,7] have made it possible 
to evaluate the details of progressive damage and fracture in composites with the use of 
established material modeling and finite element models, and considering the influence of local 
defects, through-the-thickness cracks and residual stresses. In a computational simulation, 
damage evolution quantifier such as the damage volume, exhausted damage energy, and the 
damage energy release rate (DERR) are used to quantify the structural damage tolerance at 
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different stages of degradation. Low DERR levels usually indicate that degradation takes place 
with minor resistance by the structure. Structural resistance to damage propagation is often 
dependent on structural geometry and boundary conditions as well as the applied loading and the 
state of stress. 

Some simulations [8,9] have been successful in predicting damage tolerance and failure load 
of composite structure by considering ply stresses and the corresponding stress limits for matrix 
crack growth. Thus, sensitive parameters affecting fracture toughness can be identified through 
simulation, significantly enhancing the accuracy and productivity of experiments. Energy release 
rate G and critical values of the stress intensity factors K are usually suggested as two types of 
parameters for characterizing material toughness. In this chapter, mode I interlaminar and 
intralaminar energy release rates G IC of composites are examined by computational simulation. 
Results show that computational simulation has a good predictive capability for monitoring 
damage progression in composites. 

Present research models unidirectional Graphite/Epoxy notched beam specimens using the 
CODSTRAN computational code. The effects of the orientation of the fibers with reference to 
the specimen notch direction are investigated with respect to their influences on damage and 
fracture progression characteristics. 

2.2 Methodology 

Computational simulation is implemented via the integration of three modules: (1) composite 
mechanics, (2) finite element analysis, and (3) damage progression tracking. The composite 
mechanics module [5] is designed to analyze fiber composite structures with an updated 
composite mechanics theory. Its main function is to calculate ply and composite properties of 
laminates from the fiber and matrix constituent characteristics and the composite layup. Prior to 
each finite element analysis, the ICAN module utilizes a resident data bank that contains the 
typical fiber and matrix constituent properties, computes the composite properties and 
synthesizes the laminate generalized force-displacement relations according to the composite 
layup. Additionally, ICAN detennines the composite structural response and ply stresses from 
the FEM analysis results. In simulation, the composite mechanics module is called before and 
after each finite element analysis. The finite element analysis module is capable of linear and 
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nonlinear static and dynamic analysis. Four-node anisotropic thick shell elements are usually 
used to model laminated composites [6], The finite element analysis module accepts laminate 
properties from the composite mechanics module and performs the structural analysis at each 
load increment. After structural analysis, the computed generalized node stress resultants and 
deformations are provided to the composite mechanics module. The composite mechanics 
module computes the developed ply stresses for each ply and checks for ply failure modes at 
each node. Failure criteria applied to detect ply failures are based on the maximum stress 
criterion and modified distortion energy (MDE) for combined stress effects [5], 

The overall evaluation of composite structural durability is carried out in the damage 
progression module [9] that keeps track of composite degradation for the entire structure. The 
damage progression module relies on the composite mechanics module for composite 
micromechanics, macromechanics and laminate analysis, and calls the finite element analysis 
module for global structural analysis. If excessive damage is detected, the incremental loads are 
reduced and the analysis is restarted from the previous equilibrium stage. Otherwise, if the 
increment of loads is acceptable, another finite element analysis is performed but the constitutive 
properties and the finite element mesh are updated to account for the damage and deformations 
from the last simulation. Simulation is stopped when global structural fracture is predicted. 

2.3 Microstress Level Damage Tracking 

The matrix in orthotropic composite plies is divided into two parts: regions A and B. Region 
A represents the area in which stress concentrations induced by the interaction of matrix and 
fiber do not create any effect in matrix. Region B represents the interaction zone between fiber 
and matrix. Considering the behavior of longitudinal stress <7 in , transverse stress <J l22 , in-plane 

shear stress <J n2 , out-of plane shear stress <7 m , temperature gradient AT, , and moisture M,, 
Murthy and Chamis [5] present the complete set of equations for evaluating ply microstresses in 
regions A and B. For example, ply microstresses due to <7 IU are given by: 
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in which <J mU is the matrix longitudinal stress, <7 fn is the fiber longitudinal stress, crj^ 2 is the 
matrix transverse stress in region A, and a^J 2 is the matrix transverse stress in region B. If the 

ply is subject to combined stresses, its microstresses are obtained by simply superimposing 
results of all corresponding stress components. Ply transverse fractures usually begin in region B 
due to the elevated stress levels from stress concentration. Microstress level damage tracking is 
able to quantify the type of damage in the matrix by comparison of microstresses with constituent 
stress limits. A microstress damage index is defined as a binary number with 14 bits in the 
damage progression module. The 14-bit number corresponds to the following components: (1) 
SM1A(+) longitudinal stress in region A, positive; (2) SMIA(-) longitudinal stress in region A, 
negative; (3) SM2A(+) transverse stress in region A, positive; (4) SM2A(-) transverse stress in 
region A, negative; (5) SM2B(+) transverse stress in region B, positive; (6) SM2B(-) transverse 
stress in region B, negative; (7) SM3A(+) normal stress in region A, positive; (8) SM3A(-) 
normal stress in region A, negative; (9) SM3B(+) normal stress in region B, positive; (10) 
SM3B(-) normal stress in region B, negative; (11) SM12A in-plane shear stress in region A; (12) 
SM12B in-plane shear stress in region B; (13) SM13A out-of-plane shear stress in region A; (14) 
SM23A out-of-plane shear stress in region A. When the binary bit corresponding to SM2B(+) is 
set equal to 1, it indicates that region B is fractured with the transverse failure mode. In a 
subsequent stage transverse cracking will spread to region A. 


2.4 Simulation of Braided Composite Notched Beam Specimens 


Two fiber reinforced bismaleimide prepreg coupons (X5260/G40-800) under three-point 
bending [3] are used in computational simulation. The specimen has a length of L = 27.4mm 
(1.079 in), a width of W = 6.7mm (0.264 in) and a thickness of B = 6.7mm (0.264 in). The fiber 
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volume ratio is V , =0.60 and the void volume ratio is V v =0.01. Specimens are labeled as 

‘T-type” and “JT-type” as shown in Figures la and lb, respectively. For the W -type specimen, 
the crack will grow along the fiber direction, whereas the crack will extend orthogonal to the 
fiber direction in the T -type specimen. Specimens are assumed to be dry with zero moisture 
content. The fiber and matrix properties used for computational simulation are given in Table 2.1 
and Table 2.2, and the configuration for specimens is shown in Figure 2.1a and Figure 2.1b. The 
initial crack length for T -type specimens is assumed as 2.96mm (0.117 in) and 3.08mm (0.121 
in) for W -type specimen. 

For the W -type specimen, there are 1275 nodes and 1200 Mindlin type rectangular elements. 
Its finite element model is shown in Figure 2.2b. The notch is modeled to have zero width in the 
finite element model. To represent the test setup, node at one support is restrained against 
translation to model the fixed boundary and node at the other support is constrained only 
moveable along the length of the beam. Load is applied at the top center node. Computational 
simulation indicated a damage initiation at the crack tip when loading reached 73.40N (16.5 lbs). 
The damage initiation mode was by ply transverse tensile failures in the 0° ply and the 
microstress damage indices were activated for the SM1A(+), SM2A(+), SM2B(+), SM3A(+), 
SM3A(-), SM3B(+), SM12A, SM12B, SM13A, and SM23A microfailure modes. After damage 
initiation, crack did not extend immediately. When loading was applied at 85.27N (19.17 lbs), 
crack extended again and transverse tensile failures in the 0° ply were still the only damage 
progression mode. The microstress damage indices were activated for SM1A(+), SM2A(+), 
SM2B(+), SM3A(+), SM3A(-), SM3B(+), SMI 2 A, SM12B, SM13A, and SM23A microfailure 
modes. With increased loading, transverse tensile failures continuously spread to the nodes 
surrounding the crack tip. The microstress damage indices were also activated for SM1A(+), 
SM2A(+), SM2B(+), SM3A(+), SM3A(-), SM3B(+), SM12A, SM12B, SM13A, and SM23A 
microfailure modes. Then, crack developed quickly with transverse tensile failures occurring in 
more nodes around the crack tip. The specimen continued to fracture rapidly and broke under the 
139. 54N (31.37 lbs) loading. Computational simulation indicated that transverse tensile failures 
caused fracture of the specimen. 

For the T -type specimen, we use 8 -node brick elements for constructing the finite element 
model. There are a total of 1977 nodes and 1200 elements. The finite element model is shown in 
Figure 2.2a. Only center node at one support is restrained against translation to model the fixed 
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boundary and the other nodes at the same support are constrained only to be moveable in the 
longitudinal direction of the beam. Correspondingly, center node at the other support is 
constrained only moveable in the longitudinal direction and the remaining nodes at the same 
support are constrained moveable in the longitudinal and transverse directions of the beam. 
Moreover, nodes along the crack edges are tied by duplicate node constraints to enforce uniform 
displacement. The top center nodes in the transverse direction are also considered as duplicate 
nodes to enforce the uniformly displaced loading. The damage initiation load was 37.8 IN (8.5 
lbs) by ply transverse tensile failures in the 0° ply. The damage initiation took place at the crack 
tip. The microstress damage indices were activated for SM1A(+), SM2A(+), SM2B(+), 
SM3A(+), SM3A(-), SM3B(+), SM12A, SM12B, SM13A, and SM23A microfailure modes. 
After damage initiation, transverse tensile failures spread to nodes around the crack tip. As the 
loading reached 73.04N (16.5 lbs), the crack tip developed additional damage. The microstress 
damage indices were activated for SM1A(+), SM2A(+), SM2B(+), SM3A(+), SM3A(-), 
SM3B(+), SM12A, SM12B, SM13A, and SM23A microfailure modes. Then, crack grew slowly 
with increasing loading. At 101. 69N (22.86 lbs), transverse tensile failures spread to more nodes 
around the crack tip. When the loading reached 106. 14N (23.86 lbs), the specimen broke with 
fracture propagation developing very suddenly. The microstress damage indices were also 
activated for SM1A(+), SM2A(+), SM2B(+), SM3A(+), SM3A(-), SM3B(+), SM12A, SM12B, 
SM13A, and SM23A microfailure modes. Similar to the simulation for the W -type specimens, 
transverse tensile fracture mode caused failure of T -type specimen. 

Figure 2.3 shows the relation between displacement and loading for T -type and W -type 
specimens. It indicates the failure load of T -type is smaller than that of W -type. From the 
simulation results, the failure load of IF -type is 139. 54N whereas it is 106. 14N for T -type. 
Compared with test results [3], in which the failure load for IF -type is 140N and 120N for T - 
type, it appears simulation has a good prediction for final failure load. 

Damage energy and damage volume are another two important metrics for studying damage 
progression in composites. For T -type and W -type specimens, the relationship between damage 
energy and damage volume is shown in Figure 2.4. It indicates T -type can stand more damage 
than IF -type specimen. Thus, it can be concluded that IF -type is more brittle than T -type 
specimen. 
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Figure 2.5 shows the relationship between load and damage. For W -type specimen, the 
damage increases slowly until loading reaches 108. 98N (24.5 lbs), then damage increases 
quickly representing the spread of transverse tensile failures at more nodes around crack tip. 
After reaching 132. 56N (29.8 lbs), the damage increase is relatively small until specimen breaks. 
Whereas, after loading reaches 11.12N (2.5 lbs), the damage increases uniformly with increased 
loading for T -type specimen. The plot reflects that W -type specimen is indeed much more 
brittle than T -type specimen. 

Figure 2.6 shows the relationship between the critical energy release rate Gic and crack area. 
It is observed that Gic for T -type specimen is larger than that for W -type specimen. This is 
mainly due to the fiber bridging effects that impede crack propagation in the T -type specimen. 
From the test [3], results also show that Gic for T -type specimen is larger than that for W -type 
specimen. 

2.5 Conclusions 

Experimental and computational simulation results were compared to investigate the 
progression at crack in composites. The failure patterns of the investigated fiber composite 
specimens, the critical energy release rate Gic and other available computational simulation 
results indicate that: 

1. Computational simulation, with the use of established composite mechanics and finite 
element modules, can be used to predict the failure load and fracture toughness of 
composites. 

2. Microstress level computational simulation provides the details of damage initiation, 
growth, and subsequent fracture in composites. It represents a new approach for 
investigating damage mechanisms of composites. 

3. Computational simulation is able to discern damage evolution characteristics with regard 
to fiber and crack orientations. 

4. Fracture toughness characteristics are identified by computational simulation with greater 
perspective compared to only test results. 
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5. The demonstrated procedure is flexible and applicable to all types of constituent 
materials, structural geometry, and loading. Hybrid composites, as well as laminated, 
stitched, woven, and braided composites can be simulated. 
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TABLE 2.1: X5260 Fiber Properties 


Number of fibers per end = 10000 
Fiber diameter = 0.00762 mm (0.300E' 3 in) 

Fiber Density = 1743.79 Kg/nr (0.063 lb/in ') 

Longitudinal normal modulus = 199.95 GPa (2.9E +7 psi) 

Transverse normal modulus = 19.995 GPa (2.9E +6 psi) 

Poisson's ratio ( v l2 ) = 0.300 
Poisson’s ratio ( v 23 ) = 0.45 

Shear modulus (G 12 ) = 16.548 GPa (2.40E +6 psi) 

Shear modulus (G 23 ) = 8.136 GPa (1.1 8E +6 psi) 

Longitudinal thermal expansion coefficient = 1.0E' 6 /°C (-0.55E' 6 /°F) 
Transverse thermal expansion coefficient = 1.0E' 6 /°C (0.56E' 6 /°F) 
Longitudinal heat conductivity = 301 kJ-m/hr/m 2 /°C (4.03 BTU-in/hr/in 2 /°F) 
Transverse heat conductivity = 30.1 kJ-m/hr/m 2 /°C (0.403 BTU-in/hr/in 2 /°F) 
Heat capacity = 0.712 kJ/kg/°C (0.17 BTU/lb/°F) 

Tensile strength = 2.764 GPa (400.9 ksi) 

Compressive strength = 2.419GPa (350.9 ksi) 


TABLE 2.2: G40-800 Matrix Properties 


Matrix density = 1264.94 Kg/rn 3 (0.0457 lb/in 3 ) 

Normal modulus = 4.669 GPa (720 ksi) 

Poisson’s ratio = 0.30 

Coefficient of thermal expansion = 0.77E' 4 / °C (0.428E' 4 / °F) 

Heat conductivity = 0.648 kJ-m/hr/m 2 / °C (0.868E' 2 BTU-in/hr/in 2 / °F) 
Heat capacity = 1.047 KJ/Kg/ °C (0.25 BTU/lb/ °F) 

Tensile strength = 73.706 MPa (10.69 ksi) 

Compressive strength = 55.78 MPa (8.09 ksi) 

Shear strength = 31.16 MPa (4.52 ksi) 

Allowable tensile strain = 0.02 
Allowable compressive strain = 0.05 
Allowable shear strain = 0.04 
Allowable torsional strain = 0.04 

Void conductivity = 16.8 J-m/hr/m 2 / °C (0.225 BTU-in/hr/in 2 / °F) 
Glass transition temperature = 216 °C (420°F) 
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X-Y PLANE 



Y-Z PLANE 

Figure 2.2a: Finite Element Model (Brick Element) for 7-Type Specimen 



Figure 2.2b: Finite Element Model (Shell Element) for JL-Type Specimen 
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Load vs Displacement 



Displacement (in) 


Figure 2.3: Displacement with Loading for JE-Type and E-Type Coupons 
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Figure 2.4: Damage Energy with Damage Volume for IE-Type and T-Type Coupons 
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Load vs Damage 



Load (kip) 


Figure 2.5: Loading with Damage Volume for W-Type and T-Type Coupons 


G| C vs Crack Area 



Crack Area (mm 2 ) 


Figure 2.6: Fracture Toughness with Crack Area for W-T ypc and F-Typc Coupons 
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Chapter 3 

Computational Prediction of Effective Elastic Constants 
in a Cross-Ply Laminate Under Uniaxial Loading 


The reliability and durability of composite space structures is of critical importance to assure 
safe operation in diverse environments. The main design problem is matrix microcracking at low 
temperatures that adversely affects the durability and increases the permeability of a composite 
structure. Fatigue damage accumulation due to mechanical and thermal cycling causes the 
growth of microcracks into delaminations that cause the failure of interfaces between composite 
layers. It is generally understood that thin plies microcrack less readily than thick plies, and that 
the mechanical strain required to start microcracking depends upon the stiffness of the adjacent 
"constraint" plies. Since microcracking is a function of the difference between adjacent ply 
angles, a stacking sequence chosen for a composite laminate may experience microcracking 
much more severely than suggested from laboratory tests of quasi-isotropic layups. These 
observations raise concern that the presently used methods are not adequate for the design of 
space structures that will operate reliably at cryogenic temperatures. New design methods are 
needed to take into account the effects of (1) adjacent constraint ply properties and fiber angles, 
(2) ply thickness, and (3) overall ply layup on microcracking due to thermomechanical cycling. 
This report attempts to initiate development of new composite mechanics and attendant damage 
progression/durability codes to take into account the effects of adjacent ply stiffness on 
microcrack initiation and growth. The implementation of this objective requires accounting for 
the stress concentration factors imposed by constraint plies as they affect the microcrack growth 
and evolution. Primary microcracks are typically due to transverse tensile stresses and occur as 
mode 1 dilatation cracks. The new methods account for the effect of ply thickness on microcrack 
formation and durability. The constraints imposed by adjacent plies are shown to be significant. 
The method can be used to model individual fiber-matrix interaction zones and monitor changes 
at sub-constituent levels by subzoning of each sub-slice. The approach can be used to evaluate 
the effects of adjacent ply stress concentrations on microcracking and delamination potentials. 
Thermomechanical cyclic loading efects on microcracking may also be assessed by 
incorporating the constraint ply effects. Initially, the improved method is developed for a cross- 
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ply laminate. It is expected that the method can be extended to general angle-plied laminates as 
well as braided/woven composites via the previously demonstrated macromechanics method 
based on coordinate transformations. It is hoped that the generalization of the method can be 
carried out in continuation of research. 

3.1 Effect of 0° Constraint Plies on Stiffness Degradation in Cross-Ply 
Laminates 

When cross-ply laminates are subjected to static loading in the 0° ply direction, the first 
failure generally occurs in the transverse plies. Matrix cracking is the main damage mode. 
Although matrix cracking seldom causes the final failure for laminates, it still can greatly impair 
the performance of laminates. Some tests have shown that the stiffness of laminates is reduced 
because of matrix cracking. Moreover, with matrix cracks growing, other damage modes such as 
delamination and microcracking (microcracks extend to plies adjacent to a ply that suffered 
matrix cracking) will appear. These damage mechanisms will lead to laminate failure. Thus, it is 
necessary to predict the effect of stiffness reduction in transverse plies. 

Many analytical models have been developed to evaluate stiffness response of the cracked 
laminate. Berthelot (1997) calculated stress redistribution in cracked cross-ply laminates with 
assumed longitudinal displacement fields. Abdelrahman and Nayfeh (1999) constructed 
micromechanical continuum mixture 2D and 3D models to study the stress redistribution and 
residual stiffness in orthogonally cracked laminates. McCartney (2000) predicted stress 
distribution in general symmetric laminates with uniform cracks by assuming simple 
displacement patterns that did not satisfy certain boundary conditions at the transverse crack- 
constraint ply interface. 

Pagano and Soni (1983) proposed a global/local model based on plate theory to investigate 
elastic moduli in cracked laminates. Whitney (2000, 2001) detennined effective elastic constants 
of bi-directional laminates and angle-ply laminates containing transverse cracks with Pagano ’s 
global/local model. Chattopadhyay et al. (1994, 2001) used higher-order plate theory to model 
transverse matrix cracking and local delamination. 

Leblond et al. (1996) developed 2D and 3D numerical models to calculate stiffness reduction 
in cracked cross-ply laminates. Using fracture mechanics and finite element method, Joffe et al 
(1999, 2001) analyzed stiffness response in symmetric and balanced laminates due to transverse 


NASA/CR— 2003-212351 


33 



cracking. Whitcomb (2001) employed quasi-3D and 3D finite element methods to evaluate 
material properties of cracked laminates. 

The variational method based on the principle of minimum complementary energy was first 
applied by Hashin (1985,1987) to study stiffness response of cracked cross-ply laminates under 
tensile loading. In his model, stress components only depend on x-axis direction and are constant 
across the ply thickness. Nairn et al. (1989, 1992) used the variational approach to determine 
residual stress state and predict progressive damage in cross-ply laminates. Varna and Berglund 
(1991) assumed non-uniform stress distribution across the ply thickness with variational 
methods. Praveen and Reddy (1998) applied Reddy’s layerwise theory and variational method to 
analyze stress transfer, stiffness reduction and crack opening profiles in cross-ply laminates. 
Anderssen et al. (1998) also investigated stiffness degradation in cross-ply laminates with energy 
methods, but they assumed displacement fields including crack front shape. 

In Anderssen’s model, shear stain was assumed to be zero across 0° plies. Thus, it is obvious 
that the assumption is not realistic and shear stress does not satisfy the continuity condition at the 
interface between 90° and 0° ply. In our research, we improve Anderssen’s model with assuming 
more refined displacement fields. The stress components based on new displacement fields 
satisfy boundary and continuity conditions. The geometry of laminate subjected to tensile 
loading is illustrated in Figure 3.1a. The unit cell of the cracked laminate between two adjacent 
cracks shown in Figure 3.1b is adopted for predicting effective elastic constants in the cracked 
laminate. 
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Figure 3.1a: Schematic of a [0°/90°]s Laminate with Transverse Cracks Under Uniaxial Loading 
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Figure 3.1b: One-Quarter of Unit Damaged Cell with Transverse Cracking 

3.2 Theoretical Model 
3.2.1 Lamina 1(90°) 

The displacement functions for lamina 1 in Figure 3.1b can be assumed as: 


u < ' ) (x,z) = f l (x)+<p(z)f 2 (x) 

u ( z ' } (x,z)= f 3 (x)-y (3-1) 

\ 

where fi(x),f 2 (x),f^(x) and (fiz) are unknown functions. 

The boundary conditions for displacement fields are: 

m( 1) (0,z) = 0 (3 -2a) 

m( 1) (x,0) = / 1 (x) (3 -2b) 

u®(x, 0) = 0 (3 -2c) 

From (3 -2b), we have 

0(O) = O (3-3) 
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The stress and strain fields are given as: 


e® = 


e® = 


7: 


( 1 ) 


dx 

du[' ] 

dz 

du { ' ] du[' ] 

— — + — — 
dz dx 


^.(i) _ Jd . P (i) , Ji) P (i) 

c 22 c x — *"12 c : 

crW =c W. f ( l) +c (; ) . £ -( l) 
r (l) _ Jl) . v (l) 

* xz 66 / xz 


where the elastic stiffness coefficients c n are given by 

1 

= — 

E x 

s = 

12 E z 

1 

5 22 - ~Z~ 

E_ 


’22 


’12 


"12 


"22 


'66 


5 66 


(3-4) 


(3 -4a) 
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Substituting (3-1) into (3-4), we find that 


ef = (fa)' flit) 

J\) = f 3 {*) 

7*} =</>'{z)-f 2 {x) + f'{x)-y 

h 

<4” =4’ •[/;'(*)+<*< 

h 

=^2 •|/l / W+^( Z )-/2 / W]+C 1 (l i ) '^7^ 

h 

= c 66 ■ M + fl (x) ■ f (3-5) 

M _ 

Following the symmetry and traction-free boundary conditions, stress fields are satisfied 

r!!>(o, z ) = o 

J < 7 ^ (/, z)t/z = 0 
0 
h 


\ T xz(^ z )dz = 0 (3-6) 

0 

Substituting (3-5) into (3-2) and (3-6), we know 

/,(0)=0 (3 -7a) 

/,( 0)=0 (3 -7b) 

4 * (0)=0 (3 -7c) 

/,'(<»)= 0 (3 - 7d) 

^,)-/ 2 W+|/ J 'W = 0 (3 - 7e) 

On the other hand, it is assumed that 

m( 1) (0,z) = mS 1 ) ■— (3-8) 

6 
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Thus, 


/i (•>)=“£ 


(3-9) 


where u^J is an unknown constant. 


3.2.2 Lamina 2(0°) 

The displacement field is assumed 




where u 0 2 J is an unknown constant. 


The strain-displacement and stress-strain relations are expressed as: 


r (2) _ du ,v 2> ! 

xz dz dx 
^(2) _ (2) (2) (2) (2) 

u x C 11 ~ c 12 C z 

^-(2) _ (2) (2) , (2) (2) 

u z u 12 °x ^ u 22 c z 

^?=4 2 6 V! 2) (3-n 

The expression of continuity, symmetry and boundary conditions for lamina 1 and 2 are 

u® (x, t ] ) = u i2) ( x , t x ) (3 - 12e 

) = u[ 2 \x,t x ) (3-121 

T { ^(x,t t )=T {2 J(x,t ] ) (3 -12c 

<y ? (x, t x ) = <7 (2) {x,t x ) (3 - 1 2d 

4 2) M = 0 (3-12. 

h 

Jr(; ) (/,z)t/z = 0 (3-12f 

*1 

l£>M)=0 (3-128 

<r®(x,h) = 0 (3-12h 
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Thus, we assume there exists </i{t\) and the appropriate expression for t[ 2 J(x,z) is 


■^(x,z)=T^(x,t\)- 


h(h + 2 t l ) 4 h + 2t x 3 2 

:: Z H — — Z 

f- t 2 

l 2 l 2 l 2 


(3-13) 


that is 


n 


T { ^{x,z) = y/ xz {x)- 

Wv ( X ) = C { 66 ' \<P (h )■ fl( X )+ f 3 ( X ). ] 


h{h + 2t x ) 4/7 + 2t x 3 

; ; Z H — — Z~ 


^2 ^2 


(3 -13a) 


Applying (3-10) and (3-13a) to (3-1 1) and considering (3-12), we have 
u ( ; ] {x,z) = f (x) + (j){t x )f 2 (x) -{z-t x )f' (x) 


+ 


h{h + 2t x ) (2h+t x ) 2 \ 3 

2 z 2 — z + _ r z 

ti t; t; 


h t , 


c (1) 

X= 66 


J2) 

''66 


(3-14) 


Substituting (3-10) and (3-14) into (3-11) and considering (3-12), we get 
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e[ 2) = // M + 0(*i )/ 2 ' M - ( 2 - ^ )// (*) 


+ 


(*). ■ 


hill + 2 t l ) (2/2 + t x ) 2 1 

: Z — - ~ ' 


t; 


„ Z H — — Z 

t 2 t 2 

l 2 l 2 


k\ 


£ (2) _ M 0--^ 




rl? = 4^ (0/ 2 (*)+ //(*). ' 


/z(/? + 2t x ) 4/? + 2?j 3 


„ z + — z‘ 

t~ t 2 

L 2 l 2 

(2) 


f = c<? \f( (x) + <*(», )// (*)- (z - 1, )//' (*)]+ c<? ^ 


t-, 


+ c 


WiL' 




X 


/z(/2 + 2?j ) (2 A + 1 1 ) 2 , 1 


f. 


, z +— z 
l 2 l 2 


3 


}ft , 


r ! 2) = C l ( 2 2) [// (*) + )f 2 i X )~ i Z ~ h )// (*)] + C§ ' 


A 2) 

*0 £ 


+ c 


12 




'(• 


x 


/z(/7 + 2tj ) ( 2 h + t x ) 2 1 

2 z — - " 1 


t; 


„ z + — z 

t 2 

l 2 l 2 


h 2 t , 


•L 2 ) = 4 ) k(^ 1 )/ 2 W+/ 3 / W.- 


/? (/7 + 2fj ) 4/? + 2f j 3 


„ z + — z‘ 

£ 2 t~ 

l 2 l 2 


(3-15) 


From (3-12d) and (3-12h), we obtain 



(3 -16a) 
(3 -16b) 
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Substituting (3- 16a) into (3-15) and (3- 16b), we know 

+ 2) = ui M+ (At, )/,(*)+ k,h -c,?k" W 


+c !; , #'( , i) /: W+Zj’l 1 


A (A + 2tj ) (2/z + 6 ) 2 1 


/; 


„ z + — z 

t 2 t 2 

l 2 l 2 


h 2 t , 


^.- 2) =4 2) 4^(6)-/2 w+/ 3 ui 


(•» 


/ 2(/7 + 2t|) (2/2 + 6) 2 1 3 

; 2 Z H — Z 

t 2 t 2 t 2 

l 2 l 2 l 2 


+ c n{h-z)f''(x ) 

47/ M+ 4M )// W+ c<;» ■ Ah) - 4 2 Vj" W ■ 

h 


\ r ( 2 ) ] 2 

+ _ „(2) _ ]£l2_J_ 

.(2) 


"22 


, k?r 

- ^(2) 


'22 


3.2.3 Strain Energy 

The strain energy per unit damaged cell is written as: 


t 

l 2 

h 2 t , 


t; 


(3 -17a) 


(3 -17b) 
(3 -17c) 

(3 - 17d) 

(3 -17e) 


U = )■ J J [+" • 4° + <rf ■ 4‘ 1 + ■ 4? ]* + j [<4 2) ■ ef + <4 2) • 4 2) + 4? • 4? k k 


0 0 


-P-uf’(l) 
t , J 


l 2 f , 

P = hw - p 0 
h = t\+ 1 2 


(3-18) 


in which w is the thickness of the laminate. 


Substituting (3-4) and (3-11) into (3-18), we obtain: 

6/ = J Jxj ] 2 + 2cfJ ■ 4° ■ 4 + 4 ■ k-° ] 2 + 4 ■ klr ] 2 K 


0 o 

1 a 


+ — 
2 


(Wk>- k 2) ] 


12 + 2c, ( ? • 4 2) • 4 2) + 4? • k? 1 ]’ + 4? ■ [4; 1 ] 2 }&-/>• 4 2, (/ ) 


0 


(3-19) 


NASA/CR— 2003-212351 


41 



Giving the displacement fields a small variation, then variations of displacement fields and strain 
components are defined as: 


Su^ (x, z) = 8j\ (x) + 8(p{z ) ■ f 2 (x) + <p(z)- Sf 2 (x) 
Sn ( ; ] {x,z) = 8f\ (x) + (pit, )Sf 2 (x) -(z-fi )SfP (x) 


+ . 


Su { ' ) {x,z) = Sf 3 {x)- — 


h{h + 2t i) (2h + t l ) 2 1 3 

; 2 Z H -Z 

t; t 2 t; 


h t. 


&i (2} (x,z)=$' 3 (x)+&t^ • 


z — t 


= sfi ( x ) + M z ) ■ fi ( x ) + ^( z ) ■ $fi ( x ) 

= Sfl (x) + <pf )8fi (x) -iz-t, )Sf” (x) 


+ {t x )- Sfi{x)+ df‘' {x) 

& (D = jjfe) 

Z ' fi 

§£ (2) = 

^2 

= S( P i z ) ■ fi ( x ) +</>' i z )- sfi (*) + ft* ( x ) ■ - 


h{h + 2t x ) [2h + t x ) 2 1 3 

5 z 5 z + ~ z 

t; t; t; 


nt. 


Sy {2) = (fi ) s f ( x ) + Qfl (*). ■ 


t, 


h(h + 2t x ) 4 h + 2 t x 3 - 

2 2 Z 2~ Z 

t; t; t: 


(3-20) 


From (3-3), (3-7a), (3 -7b), (3-7e) and (3-9), it is obvious that 

< y ,( o)=o 

< r ! ( o)=o 

( o )=&2 

h 

5(p{ 0) = 0 (3-21) 
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The variation of energy U is given as: 


L h 


su =] dx \ {4 ■ 4 ■ 4° + 4 ■ [4 ■ + 4 ■ &?]+ 4 ■ 4 ■ 4° + 4 ■ 4 


0 0 


f , f k ( ? -4 2) •*?' +4? •[•£•?) •*?' +ff -&f ]+ 4 > f‘ 2) -&? 1 
! 1,1+ 

-P.<S7,' 2 »(/) 


•fife 


Substituting (3-14) into (3-16) and integrating it over z-direction, we have 
SU = SU (1) + SU (2) 


where 


<St/ (1) = J 
0 



h h h 

+ 1 6 (x)J S(p{z)dz + I 7 (x)J <p(z)S<p(z)dz - 1 8 (x)J <p (z )S<p(z)dz 

0 0 0 


fife 


Ij (x) 4^1 




fi ( x ) ■ J \<t> ( z )} 2 dz + 4 ■ 4| ) ■ // (*) 




1 3 M = 4 ■ // M ■ J <44 + 4 • fi ( x ) • {4 ( z 4 + 4 • 

0 0 




! 5 W 

i 6 W 

I 7 (x) 

! 8 W 



■ dyii}dz 

(3-22) 

(3-23) 

(3 -23a) 


NASA/CR— 2003-212351 


43 



and 


SU (2) = | 


1 9 W + I.O ( X )<f 2 M + 1 11 {xffll (*) + I 12 (*Wl (*) +I 13 ( X )C / (*) 

+ I 14 (x)<$/ ( (2) 


*-Oz 


dx 

(3 -23b) 


I 9 (*) = 4?* 2 • fi (*) + 4? ] t 2 • ) ■ / 2 (x) + cj? -y- • <P‘ {t x ) ■ fl (x)-c} ? y ■ f" (x) 


At l 


+ c n , T^-/3 // W + c i ( 2 2) -Wo? 


1 10 (-4 : 


,( 2 ) 

''66 


12 

2A 

15 


-k(6)] 2 - 


66 •//(*) 




111 (*) = 4l4 ■ ) ■ fi (x) + cff ■ <p' (fj ) ■ // (x) + 4 2 4 • f- {t x ) ■ fl (x) 


+ C 


+ C 


(2) £t\ 


■ k (6 )]“ ■ // (*) + 4?' Ah ) ■ fl (x) - cl? k ^ ) ■ fl' (x) 


105 


(2) 


12 




2.3 


■ ) ■ fi" i x )~ 4? 44“ ' ^ (0 ' // (*) + cl? ■ <p‘ (fj ) ■ f?' (x) 


30 


105 


+ c| 2 2) ■ (/>(t x ) ■ uJJ + cl? . u (J 


2/1 2 ? 


2/l 2 ^ 


T i 2 (*) = 4? -yy ■ fi ih ) ■ / 2 4) + 4? -yy ■ fl (x) 

J i3 04 = 4? y L ■ fi 04“ 4? 1 y ■ fl (x) + 4? y- ■ 0(f, ) ■ fl (x) - 4 2) y ■ </>{h ) ■ fi 04 

I'M 


2.3 


■ 4? ^ • // to + 4? ^ •/, (x) - cj? | • u$ 


+ c 


12 


■ u 


(2) 

Or 


1 14 (x) = 4? • fi 04 + 4? ■ ^Oi )/ 2 4) + c i ( 2 2) y- ■ fi ' Oi ) ■ fl to - - C l? ■ y f? (x) 


+ C 


l (2) Tk-/3 // W + 4 2) — -«0? 


12 
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Integrating (3-23a) and (3-23b) again by parts and considering (3-21), it yields 


SU = J {[- or, ■ f(' (x) - a 2 ■ ./; (x) - co x ■ f' (x) - co 2 • / 3 (3) (x)]#i {x)}dx 
0 
L 

+ J 1“ «2 • fl ( X ) + «3 ' / 2 (*) - «4 ' f- (*) + a 5 ■ fl (*) - a i ■ /i 3) (*)]#2 ( X )K 


+ 


J- 


■ fl {x) + (0 6 ■ fi [3) {x) + a 9 ■ f 2 (x) + ffl 10 • /; (3) (x) + ■ f 3 (x) 

-6\- f^{x)+CO u -f^(x) 


Sf 2 {x)\a 


+ Jte +Pi -</>{ z )-p3 {z^<p{z)}dz 

0 

+ [«,/,' (/)+«,//(')+ ®,/, (/)+ ©,// w + «,«£’ - %, ( i ) 

+ [Iff, + )/| (/) + (ft 4 + j )_/~2 (/) + (/) + (ff 7 + )_/j (/) + — (O^ 

+ k/i (o) + a>, 0 X," (0) + £W, , X, 1 1 1 (0)]<Sf , ( ,'i +k/,(0+‘» 1 /2(Z)+‘»5"l?]*S? 


where 



% 2 (/) 

(3-24) 
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a, 


= c[f A ■ (p 2 (?,)+ c, (2) ^ A (O' (pit,)- eg ^ ■ (p 2 [t, ) 


At; 


„( 2 ) 


2,3 


15? 


6t , 


c? = A 2 ^ 

m 12 ‘'ll 


2^ 

15?, 




,( 2 ) 


2,3 


2/1 ?_ 
105?, 


•^i) 


2A t_ 

105?, 


?y, c,2 

®2=q ( ?— "^)- 
11 12 2 

«3 = C!? 

co A = c{ 2 2) ■ (pit, ) + c£> ■ (p j (?, ) + cf? tjL ■ ^(?, ) - eg ^ ■ ^(?, ) 

12 ?, 6 ?, 


[ i 


.( 2 ) 

'22 


/ 


ft) = r^t~ 

U/( L,J 1 2 


A_I 

12 2 


.(i) 


= eg ■ ^(?i )- — f (piz)dz + 4 2) ■ <p' (?, ) 


^1 0 

( 1 ) ?, ( 2 ) 2/l 2 ? 2 

®8=egy + cg ) -^- 

?y 9 = eg — 


15 


? 


?y, 0 c,[ 


I? ^ ) + Cn ■ <t>' ih )~ eg y ' )" c® ^ ■ (p' (?, ) 

v 3 


2.3 


=cgA + c (2)Al_ c |2)Ai 

11 11 3 11 105 11 15 

?y , 2 = (pit l ) + — ■ (pit,)-^- Pit x ) + ^z- ■ (p ! it,) 
?, 6 ?, 12 


and 


L (l) L (l) L 

A = eg } [/; (x) ■ fi (x)J/x + ^ J |/. (x) ■ / 3 (x)J& - C -J- 1 1 f 2 (x) ■ // (x)}/; 

0 'l 0 M 0 


■<p'itl)-<pitl) 

(3 -24a) 


(3 -24b) 


NASA/CR— 2003-212351 


46 




0 


Pi =4 1 6 ) J[/' 2 W] dx 

0 

According to the theorem of minimum potential energy, there is 

SU = 0 


So, for all admissible displacement, there are 

«1 ■ fx U) + oc 2 ■ f” (x)+&>, ■ f' {x)+(D 2 - / 3 (3) (x) = 0 
«2 -fi(x)-a 3 -f 2 (x) + a 4 - f" (x)-a 5 ■ f' (x) + a 7 • / 3 (3) (x) = 0 
' fx (•*•) ’ f\ i X ) + ^9 4 f 2 (•*■ ) ®10 ' fl ( X ) ^9 ' f 2 (•*') 

-® 8 -//W+^n •/ 3 (4) W=° 

P\ +Pl ■ fa)- Pi = 

A ' ( 0 + <*2/2 ( 0 + ®i / 3 ( 0 + ^2/3" ( 0 + ® 3 «o? - f 5 = 0 
(cr 2 + CTjo )y i (/) + [ot 4 + cx n )f 2 (/) + oc b f 2 (/) + ( oc 7 + cc n )y 3 (/) + cc % u\ )z — co n 
cojf (/) ■ - (oJ 2 (/) ■ - fflg// (/) + t^yf (/) + co n f!f ] (/) = 0 
^// / (o) + ® 10 y2 // (o)+®n/3 (3) (o) = o 
® 3 yi(0+®4/2(0+®s“o? =0 

Solving (3-29), we obtain 

<pf) = \ sinh(^z) + B $ cosh (k $ z) - C , 

where 



From (3-3) and (3-7c), there are 

4 =0 

B o=C, 


(3 -24c) 


(3-25) 


(3-26) 

(3-27) 

(3-28) 
(3-29) 
(3-30) 
P = 0 (3-31) 

(3-32) 
(3-33) 

(3-34) 


(3-35) 


(3 -35a) 
(3 -35b) 


(3-36) 


NASA/CR— 2003-212351 


47 



Thus, 


0(*) = C,[cosh(A,z)-l] 


From (3- 17c), we know 


fi (*) = CJi (x) + ) ■ fi (x) + / 4 / 3 (x) 


C: 


,(i) 

"12 


3 c (2) t 

*"12 *2 

c (1) 

Z' — ii 

c (2) tt 
*"12 ^r2 


After substituting (3-38) into (3-26), (3-27) and (3-28), we have 

S/Z (*) + S/ 2 " (*) + K zfi (*) = 0 

*■ 4 / 1 ' " (*) - a zfi W + S fi (x) + K Jz (*) = 0 

K ifl (*) + s/ (3) (*) + s// (*) + So/i 3) (*) + S1/3 (*) = 0 


where 


S =«, +^ 2 / 

S = a 2 +oj 2 C 3 -(p^) 

k 2 = OJ ] + 00 4 

S =a 2 +a 7 C, 
tc 5 =a 4 + a 7 C 3 ■ 0(f t ) 

S — o^iC 4 ■ 

S — <W n ^ 3 ^ 4 - <^8^3 

S =«6+«11^3 

S —0^9 +<^ 11 / 3/4 ^(^1 ) “ ^8^3 '^(h) 

50 = ®L0 ^1 ! £3 ' ) 

51 =C0 9 +C0 n C 2 -(O s C 4 

Then, there are two new equations 

Yifi" (*) + Yifi ( x ) + Yzfi (*) = 0 

yJ " (*)- Y 5 f (4] {x)+ yj 2 (x)- 7 7 /, (4) (x) = 0 

in which 


(3-37) 


(3-38) 


(3-39) 

(3-40) 

(3-41) 


(3 - 41a) 


(3-42) 
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y i = ^1^6 ~ k ^ k a 

Yl = **3«3 

r 3 = *2 *6 - *3 *5 

y 4 =*’i*‘ii 

r 5 = 

7 6 ~ K 2 K \ \ ~*7S 

7i=^0 (3 -42a) 

From (3-42), we obtain 

^,/ 2 (4) (*) + ^ 2 /)> (2) (*) + Vji (x) = 0 (3-43) 

where 

V'i = ^3 r 5 - ^7 

Wi = r 2 r 5 + rir 6 - r 3 7 4 

¥i =7 i7a (3 -43a) 

Solving (3-43) and considering (3-6) and (3-12) as well as symmetry, we have: 
fi {x) = d l sinh(cKr) + d 2 sinh(yfir) (3-44) 

for 


^<0 

¥i 

vl~*YiYi *0 


where d\ and <7 2 are constants and 



Substituting (3-44) into (3-42) and considering (3-7), we obtain 

/) (jc) = d 2 x - + d x ■ sinh(cpc) - ^ 2 + d 2 ■ sinh(^r) 

7i« YiP~ 


(3 -44a) 


(3 -44b) 


(3-45) 
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where ck is a constant. 


Then, we can know d?, and fi(x) from (3-7), (3-9) and (3-41) 

d i 


<f 3 = 


A i-i 2 i i Yi K * a i K ^ a 




Yx K i 


Yx K i 


K n 


K 


+ 


7 J 
3 A 


Yi \ ft l I MsP_ } %P~ k*P ^xoP 


YxP 7x 


Yx K i 


Yx K i 


Kn 


K 


d-, — 


7 J 


2 “Or 

Kn 


f 3 (*) = 


d x [cosh(ew) - 1] 


+ 


Yi^i +Yd^0(_ +jw*_ jW -YxW' -YxKw<*' 

Yx K xx® 

Mi + -MiE !Mi£ zMi£ ^[coshW-ii-H 


Yi*uP 


,W 

0 z 


(3-46) 


With Newton-Raphson method, we can solve a and j5 from (3-35a), (3-44a) and (3-44b). Then, 
substituting (3-44), (3-45) and (3-46) into (3-30), (3-31), (3-32), (3-33) and (3-34), it is obvious 
that (3-33) is satisfied automatically. Thus, solving equations (3-30), (3-31), (3-32) and (3-34), 
we obtain parameters d\, di, u^} and u\f} . 


3.3 Effective Moduli 


All average strain and stress components are calculated with: 

h l 


g (1) =- i -\dz\g [l] dx 

M ^ n a 


1 0 0 
h L 


g (2) =^-\dz\g [2) dx 


M 0 


(3-47) 


Thus, the effective elastic modulus E x is: 

^(0 

F W - 


C x 

M2) 

7(2) _ °Y_ 
'* £ & 


(3-48) 
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3.4 Conclusions 


A new analytical method with spatially consistent continuity conditions has been developed 
for the evaluation of stiffness degradation of cross-ply composites with transverse matrix cracks. 
The developed method is able to quantify the correct residual stiffness of transversely cracked 
plies, taking into account the effects of adjacent constraint ply properties. The development of 
this method is an essential step to improve accuracy of progressive fracture simulations of 
composite space structures. The method can be extended to general angle plied composites via 
macromechanics and coordinate transfonnations. 
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Chapter 4 

Calibration and Validation of Software Program IC AN/JAVA 


4.1 Objective 

The objective of this chapter is to present the comparison of results using the old ICAN (4) 
version and the new ICAN/JAVA (8) version. The research also presents the comparison of 
results using the Multi Factor Interaction model not available in the old ICAN. 

4.2 Background 

The most cost effective way to analyze/design composite structures is through the use of 
computer codes. Over the last two decades the research in composite micromechanics and 
macromechanics, which includes the effects of temperature and moisture, has resulted in the 
development of several computer codes for composite mechanics and structural analysis. 

The need for a multilevel analysis to design structural components made of fiber composites 
is necessary for better design. The Multilayered analysis consists of (1) Micromechanical 
theories for the thennoelastic properties and the stress level limit of the single ply as function of 
constituent material properties and the particular fabrication process, (2) the combined stress- 
strength criteria for the single ply, and (3) Multilayered composite structural response and 
analysis where the interply layered effects are taken into account. This code is identified as 
MFCA(l). 

Intraply hybrid composites are a logical sequel to conventional and interply hybrid 
composites. Another code INHYD (2) integrated the mechanical behavior of Intraply hybrids 
with those for hygrothermal effects for predicting hygral, thennal, and mechanical properties of 
Intraply hybrid composites for their design. 

ICAN is the synergistic combination of the micro mechanical design of INHYD and the 
laminate analysis of MFCA with several significant enhancements. It was primarily designed to 
analyze the hygro-thermo-electro-mechanical properties and response of fiber- or particulate- 
reinforced, resin-matrix or metal-matrix layered composites, given the local membrane loads and 
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bending moments. Three types of layers are recognized by the program: (1) the standard 
composite system that consists entirely of a primary composite made of one type of fiber and 
matrix; (2) the intraply hybrid composite system that consists of a primary composite and a 
secondary composite arranged in a prescribed manner within a layer (For purposes of 
identification, the primary composite in the hybrid is the one that constitutes the largest volume 
ratio.); and (3) the interply layer that consists of the matrix. In addition, ICAN recognizes 
moisture, temperature, and electrical gradients through the thickness. However, within each ply 
(or slice) the temperature, moisture, or voltage is assumed to be constant. 

During the last six years, since its intial release ICAN has undergone many changes to fix the 
bugs that have surfaced and also improvement to its micromechanical equations. ICAN/JAVA is 
the result of these enhancements which had made it more manageable and user-friendly. 

4.3 Features and Enhancements of ICAN/JAVA 

1 . The output can now be tailored to specific needs by choosing the appropriate options. 

2. Several modules have been added to perform durability/fatigue type analyses for thermal 
as well as mechanical cyclic loads. The code can currently assess degradation due to 
mechanical and thermal cyclic loads with or without a defect. 

3. The laminate configuration is not restricted to only plies but can be sliced and subsliced 
for a closer look at what goes on in the ply layer. 

4. The laminate configuration can also be used to incorporate any defects in the layup. 

5. The thermal loads, hygral loads and electrical loads can now be input as constant, linear, 
parabolic, hyperbolic or user defined across the ply layup. 

6. The damping analysis has been incorporated in the new version. 

7. Details regarding the reaction of the metal and matrix composite can also be input in the 
new version. 

8. Details regarding the impact of a hard spherical projectile crashing into the composite can 
also be input in this version. 

9. The ability to change empirical constants like the experimental correlation factor for 
combined stress value criteria, and for onset of delamination criteria value is also 
provided. 
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4.4 Theory 


In ICAN/JAVA the initiation of damage is detected using modified distortion energy (MDE) 
failure criteria (33) The MDE failure criteria is a variation of the Tsai-Hill theory (34,35) that 
states that failure is initiated when the following inequality is violated. 



Where K L i 2 is the longitudinal-transverse directional elastic properties interaction factor, 
given in terms of the lamina elastic constants by 


(1 + 4v Ll2 

^L\3 )E L21 + (j_ ^ L23 )e l ii 

IE L uE L22 {2 + V l \2 +^113 )(2 + Fl 21 + ' / I 23) 


a L n and, ol 22 are the lamina stresses along and perpendicular to the fibers and au 2 s is the shear 
stress and F is the failure. 

Sui and S L 22 , are the lamina strengths along and perpendicular to the fibers and S L i 2 s is the 
shear strength. If the prevailing stress in the fiber direction is tensile (aui a >0) then the strength 
Sui is the tensile longitudanal strength Slut and if the prevailing stress in the fiber direction is 
compressive (ou i a <0) then the strength S L n is the compressive longitudanal strength S L nc- The 
same is applied for the transverse stresses. 

E L ii and E L 22 are the lamina moduli along and perpendicular to the fibres Vli 2 , Vlb and Vl 23 are 
the lamina poisson’s ratios 

Failure is deemed to take place when F=0 and failure modes can be either fiber or matrix 
failure depending upon the dominanat terms in the MDE equation. Fiber failure is assigned when 
one of the following equations is satisfied 


NASA/CR— 2003-212351 


55 



(4.3) 


LWa 


V ^LUa J 
or 

& LWa 

V ^LUa J 


> 


> 


^ Lllb 

V ^ LI lb J 


^ L\ Is 

V ^Ll l.s J 


The stress limits Slut, Sluc and Sli 2 s in the MDE failure criteria are calculated by ICAN 
based on the constituent fiber and matrix strength and micromechanics equations. The lamina 
elastic properties used in the directional interaction factor K L i 2 are also computed by ICAN from 
the fiber and matrix elastic properties via micromechanics equations (4). 

The simulation of complex material behavior resulting from the interaction of several factors 
such as temperature, fatigue, time dependence etc has been mainly performed by factor-specific 
representations like creep analysis, fatigue analysis etc. Suppose we assume material behavior to 
be a continum represented by some surface (Reference 12). Initially a primitive fonn of MFIM 
representation for describing complex composite behavior in polymer matrix composites (36). 
This was extended to metal matrix composites (37) and is continuing. 

A multifactor interaction model is briefly described to represent complex point material 
behavior in a single equation. The model is of product form in order to represent coupled 
interactions and to be computationally effective. The model describes a continuum or surface in 
space that represents the complex material behavior in terms of various factors that affect a 
specified material behavior as shown in equation below (38). 


M t 


M 


po 


f t* j 1 \ 

1 gw 1 

T -T 

\ A gw 1 o J 


at 

Sftf j 
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G ' M 

SfNjM j 


a T N T 
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v SfNjj. j 


(4.4) 
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<4 J 
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of J 


where M=property, T=temperature, S=strength, o=stress, N=number of cycles, t=time, co=load 
frequency, E e =erosion effect, C e =corrosion effect, C c =chemical/metallurgieal effect 
Subscripts: gw=wet glass temperature, o=reference condition-assumed nominal at ambient 
conditions, f=flnal condition, M=mechanical load, T=thermal cyclic load, M,n,etc are exponents 
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for the mterial that property effect which describe respective behavior paths from the reference to 
the final value. 

ICAN/JAVA has integrated the following factors in its code: thermal, stress, stress rate, 
temperature rate, reaction, mechanical cycle fatigue, thennal cycle fatigue, time at stress, 
moisture, porosity, hygrothermal effect and frequency effects. 

Reference value, final value, increasing exponent, decreasing exponent, and other exponent are 
entered in the material bank of ICAN/JAVA. 

The equation of the multiplicative factor is 


final _ value - factor _ value 
final _value - reference _value 


~i exp orient 


(4.5) 


where the program chooses the exponent in the computation from the three exponents 
(increasing exponent, decreasing exponent and other exponent) as appropriate. For the 
hygrothermal effect the equation of the multiplicative factor is slightly modified as below 
(Reference 14): 


T -T 

gw 


L gd 


exp onent 


(4.6) 


where T 0 is the room temperature with a default value of 70 degrees Fahrenheit, T g d is the dry 
glass transition temperature, entered as the final value, T gw is the wet glass transition 
temperature. T gw is not entered but is computed from T g d and the moisture at the current layer as 
below: 

T gw =T gd (l-0.1 M f +0.005 M/) (4.7) 

where M f is the moisture fraction, T is the computed temperature value of the current 
ply/layer/slice. 
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4.5 Model Description 


4.5.1 Durability Analysis 

The laminate consists of 16 plies as shown in Figure 4.1 with a ply thickness of 0.005 inch. 


0 DEGREE PLY 

16 

45 DEGREE PLY 

15 

-45 DEGREE PLY 

14 

90 DEGREE PLY 

13 

0 DEGREE PLY 

12 

45 DEGREE PLY 

11 

-45 DEGREE PLY 

10 

90 DEGREE PLY 

9 

90 DEGREE PLY 

8 

-45 DEGREE PLY 

7 

45 DEGREE PLY 

6 

0 DEGREE PLY 

5 

90 DEGREE PLY 

4 

-45 DEGREE PLY 

3 

45 DEGREE PLY 

2 

0 DEGREE PLY 

1 


Figure 4.1: Composite Ply Lay-up for Durability Analysis 

The material is made of T300/IMHS composite the properties of which are shown in Chapter 
5.10. The fiber volume ratio for the primary material is 0.55 while for the secondary material is 
0.54. The cure temperature (Tcu) is 160F while the use temperature (Tu) is 350F. The moisture 
content is 0.5%. 

The laminate is subjected to in-plane membrane loads of 100 pounds per inch and 50 pounds 
per inch in the x and y directions. Two different cases are run and the results compared for 
durability. In addition to the static loads in the first case the tension-tension cyclic load of 200 
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pounds per inch is input. The cyclic degradation coefficient is taken as 0.01 and the number of 
cycles is 10 million. 

In the second case the tension-tension and the bending cyclic load of 200 pounds inch per 
inch is input with the same cyclic degradation coefficient and number of cycles as shown in 
Table 4.2. 


Table 4.1: Durability Analysis Loads for the Test Cases 


Mnemonic 

Test Case 1 

Test Case 2 

Nxx 

100 

100 

Nyy 

50 

50 

CNxx 

200 

0 

CMxx 

0 

200 

Cyclic Degradation Coefficient 

0.01 

0.01 

Number of Cycles 

l.E+07 

l.E+07 


This problem is compared with ICAN, ICAN/JAVA (subslicing=0) and ICAN/JAVA 
(subslicing = 9). The BOOLEAN MICRO is set true because the ply properties were not found in 
the online databank for ICAN/JAVA. Hence to compare correctly both the MICRO keywords are 
set to true in the input files. 


4.5.2 MFIM Analysis 

The laminate consists of 4 plies as shown in Figure 4.2 with a thickness of 0.01 inch for zero 
degree plies and 0.005 inch for ninety degree plies. 


0 DEGREE PLY 

4 

90 DEGREE PLY 

3 

90 DEGREE PLY 

2 

0 DEGREE PLY 

1 


Figure 4.2: Composite Ply Layup for MFIM Analysis 


It is made of AS— /IMLS and SGLA/HMHS (Chapter 5.10) composite materials with the first 
applied to the zero degree plies and the second to the ninty degree plies. The primary fiber 
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volume ratio for AS— /IMLS and SGLA/HMHS is 0.55 and the secondary fiber volume ratio is 
0.57. The cure and use temperature is 70F. 

The laminate is subjected to in-plane static load of 1000. pounds per inch in the x-direction 
and a combination of cyclic loads shown in Table 4.1. 


Table 4.2: Loads Applied for MFIM Test 


Mnemonic 

Upper limit of 
cyclic load 

Lower limit of 
cyclic load 

Number of 
cycles 

Cyclic degradation 
coefficient (Betal) 

CNXX 

200. 

100. 

100. 

0.1 

CNYY 

-50. 

-100. 

10. 

0.1 

CNXY 

20. 

10. 

100. 

0.2 

CMXX 

10. 

5. 

10. 

0.01 

CMYY 

4. 

2. 

1000. 

0.15 

CMXY 

2. 

1. 

100. 

0.1 

NXX 

1000 

- 

- 

- 


The material properties are further changed with the help of the MFIM using the following 
format as shown in Table 4.3. 


Table 4.3: MFIM Input Options as Shown in ICAN/JAVA 


Computed Factor 

Reference 

Value 

Final Value 

Increasing 

Exponent 

Decreasing 

Exponent 

Stress (psi) 



0.5 

0.5 

Stress rate (psi/sec) 



0.5 

0.5 

Temperature rate (F/sec) 

73 

73 

0.34 

0.34 

Mechanical Cycle Fatigue (hertz) 



0.5 

0.5 

Time at stress (sec) 



0.5 

0.5 


The material properties that were changed were the Poisson’s ratio, the elastic moduli and the 
strengths. The new composite was named AS-1/IML1 and SGL1/HMHS and the same 
simulation carried out and the results compared. 


4.6 Results and Discussion 

The durability results are compared for the different loads. The graphs are plotted between 
the safety factors and the ply layup. A brief explanation of these safety factors is as follows. Mil 
is a safety factor defined as below: 
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Mil = 1 — 1 


whereRl 1 = — (4 ’ 8) 

c nr c 

° LUT u ' ‘-’LUC 

where R1 1 is defined as the ratio of the stress in 1 1 direction to the ultimate stress in tension or 
compression. 

M22 is defined below: 


M22 = \-R22 

whereR22 = ^ (4 ’ 9) 

9 nr 9 

^ L22T U ' lJ L22C 

where R22 is defined as the ratio of the stress in 22 direction to the ultimate stress in tension or 
compression. 

M 12 is a safety factor defined as below: 


M12 = \- R\2 
whereR\2 = ^ 


12 


(4.10) 


’L\2S 


where R12 is defined as the ratio of the stress in 12 direction to the ultimate shear stress. 


4.6.1 Cyclic In Plane Membrane Load 

Cyclic in plane membrane loads were input and the results computed from ICAN (Case 1), 
ICAN/JAVA (sub-slicing=0) (Case 2) and ICAN/JAVA (sub-slicing=9) (Case 3) 

The factor of safety for the transverse, longitudinal and shear directions were calculated and 
enumerated in Tables 4. 4-4. 6 and the variations plotted with respect to ply layup as shown in 
Figures 4. 3-4. 5. 
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Table 4.4: Factor of Safety in the Longitudanal Direction for Cyclic in Plane Membrane Load 


PLY 

ANGLE 

Mil 

ICAN 

ICAN/JAVA 
(Without sub-slicing) 

ICAN/JAVA 
(With sub-slicing) 

1 

0 

0.934 

0.982 

0.982 

2 

45 

0.976 

0.982 

0.982 

3 

-45 

0.976 

0.982 

0.982 

4 

90 

0.930 

0.967 

0.967 

5 

0 

0.934 

0.982 

0.982 

6 

45 

0.976 

0.982 

0.982 

7 

-45 

0.976 

0.982 

0.982 

8 

90 

0.930 

0.967 

0.967 

9 

90 

0.930 

0.967 

0.967 

10 

-45 

0.976 

0.982 

0.982 

11 

45 

0.976 

0.982 

0.982 

12 

0 

0.934 

0.982 

0.982 

13 

90 

0.930 

0.967 

0.967 

14 

-45 

0.976 

0.982 

0.982 

15 

45 

0.976 

0.982 

0.982 

16 

0 

0.934 

0.982 

0.982 



PLY-LAYUP 


Figure 4.3: Graph of Ml 1 for Cyclic in Plane Membrane Load 
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Table 4.5: Factor of Safety in the Transverse Direction for Cyclic in Plane Membrane Load 


PLY 

ANGLE 

M2 2 

ICAN 

ICAN/JAVA 
(Without sub-slicing) 

ICAN/JAVA 
(With sub-slicing) 

1 

0 

0.666 

0.790 

0.790 

2 

45 

0.640 

0.774 

0.774 

3 

-45 

0.640 

0.774 

0.774 

4 

90 

0.614 

0.757 

0.757 

5 

0 

0.666 

0.790 

0.790 

6 

45 

0.640 

0.774 

0.774 

7 

-45 

0.640 

0.774 

0.774 

8 

90 

0.614 

0.757 

0.757 

9 

90 

0.614 

0.757 

0.757 

10 

-45 

0.640 

0.774 

0.774 

11 

45 

0.640 

0.774 

0.774 

12 

0 

0.666 

0.790 

0.790 

13 

90 

0.614 

0.757 

0.757 

14 

-45 

0.640 

0.774 

0.774 

15 

45 

0.640 

0.774 

0.774 

16 

0 

0.666 

0.790 

0.790 
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Figure 4.4: Graph of M22 for Cyclic in Plane Membrane Load 
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Table 4.6: Factor of Safety in the Shear Direction for Cyclic in Plane Membrane Load 


PLY 

ANGLE 

M12 

ICAN 

ICAN/JAVA 
(Without sub-slicing) 

ICAN/JAVA 
(With sub-slicing) 

1 

0 

1.000 

1.000 

1.000 

2 

45 

0.966 

0.975 

0.975 

3 

-45 

0.966 

0.975 

0.975 

4 

90 

1.000 

1.000 

1.000 

5 

0 

1.000 

1.000 

1.000 

6 

45 

0.966 

0.975 

0.975 

7 

-45 

0.966 

0.975 

0.975 

8 

90 

1.000 

1.000 

1.000 

9 

90 

1.000 

1.000 

1.000 

10 

-45 

0.966 

0.975 

0.975 

11 

45 

0.966 

0.975 

0.975 

12 

0 

1.000 

1.000 

1.000 

13 

90 

1.000 

1.000 

1.000 

14 

-45 

0.966 

0.975 

0.975 

15 

45 

0.966 

0.975 

0.975 

16 

0 

1.000 

1.000 

1.000 



Figure 4.5: Graph of M12 for Cyclic in Plane Membrane Load 
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From the figures we observe that plies 4 and 13 are the first to fail in the transverse and 
longitudinal directions, which is consistent with both ICAN/JAVA and ICAN. But in case of 
shear the ±45 plies are the first to fail and the 0 and 90 plies do not take the shear loads at all. 
This is consistent with the loads applied. Also we observe that there is a marked difference 
between the factors of safety for ICAN/JAVA and ICAN. This can be due to the sub-slicing used 
in ICAN/JAVA, which enables the code to predict the stresses in more detail and accuracy. 
Hence the factors of safety are usually higher or equal in ICAN/JAVA than in ICAN. 


4.6.2 Cyclic In Plane Bending Load 

Cyclic in plane bending loads were input and the results computed from ICAN (Case 1), 
ICAN/JAVA (sub-slicing=0) (Case 2) and ICAN/JAVA (sub-slicing=9) (Case 3) 

The factor of safety for the transverse, longitudinal and shear directions were calculated and 
enumerated in Tables 4. 7-4. 9 and the variation plotted with respect to ply layup as shown in 
Figures 4. 6-4. 8. 


Table 4.7: Factor of Safety in the Longitudanal Direction for in Plane Bending Load 


PLY 

ANGLE 

M2 2 

ICAN 

ICAN/JAVA 
(Without sub-slicing) 

ICAN/JAVA 
(With sub-slicing) 
Slice 1 

ICAN/JAVA 
(With sub-slicing) 
Slice 125 

1 

0 

-4.33 

-1.594 

-1.766 

-1.423 

2 

45 

-0.15 

0.427 

0.384 

0.470 

3 

-45 

-0.32 

0.341 

0.283 

0.399 

4 

90 

0.032 

0.541 

0.488 

0.594 

5 

0 

-1.49 

-0.221 

-0.392 

-0.049 

6 

45 

0.543 

0.768 

0.726 

0.811 

7 

-45 

0.622 

0.807 

0.749 

0.865 

8 

90 

0.927 

0.965 

0.912 

0.982 

9 

90 

0.811 

0.911 

0.981 

0.840 

10 

-45 

0.760 

0.888 

0.931 

0.844 

11 

45 

0.701 

0.859 

0.891 

0.827 

12 

0 

-0.84 

0.117 

0.246 

-0.011 

13 

90 

-0.39 

0.344 

0.415 

0.273 

14 

-45 

0.055 

0.538 

0.582 

0.495 

15 

45 

0.185 

0.603 

0.635 

0.571 

16 

0 

-2.96 

-0.912 

-0.784 

-1.041 
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Figure 4.6: Graph of Ml 1 for Cyclic in Plane Bending Load 


Table 4.8: Factor Of Safety in the Transverse Direction for in Plane Bending Load 


PLY 

ANGLE 

M2 2 

ICAN 

ICAN/JAVA 
(Without sub- 
slicing) 

ICAN/JAVA 
(With sub-slicing) 
Slice 1 

ICAN/JAVA 
(With sub-slicing) 
Slice 125 

1 

0 

0.352 

0.493 

0.473 

0.513 

2 

45 

0.661 

0.790 

0.790 

0.790 

3 

-45 

0.661 

0.790 

0.790 

0.790 

4 

90 

0.656 

0.785 

0.738 

0.790 

5 

0 

0.519 

0.651 

0.631 

0.671 

6 

45 

0.661 

0.790 

0.790 

0.790 

7 

-45 

0.661 

0.790 

0.790 

0.790 

8 

90 

0.656 

0.790 

0.790 

0.790 

9 

90 

0.515 

0.679 

0.789 

0.569 

10 

-45 

0.496 

0.666 

0.707 

0.626 

11 

45 

0.336 

0.543 

0.592 

0.493 

12 

0 

0.666 

0.790 

0.790 

0.790 

13 

90 

- 0.61 

- 0.20 

- 0.092 

- 0.312 

14 

-45 

0.058 

0.338 

0.379 

0.297 

15 

45 

- 0.18 

0.148 

0.197 

0.098 

16 

0 

0.666 

0.790 

0.790 

0.790 
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Figure 4.7: Graph of M22 for Cyclic in Plane Bending Load 


Table 4.9: Factor of Safety in the Shear Direction for in Plane Bending Load 


PLY 

ANGLE 

M2 2 

ICAN 

ICAN/JAVA 
(without subslicing) 

ICAN/JAVA 
(with subslicing) 
Slice 1 

ICAN/JAVA 
(with subslicing) 
Slice 125 

1 

0 

0.900 

0.909 

0.903 

0.915 

2 

45 

-0.37 

-0.226 

-0.320 

-0.132 

3 

-45 

-0.16 

-0.037 

-0.131 

0.057 

4 

90 

0.940 

0.946 

0.939 

0.952 

5 

0 

0.953 

0.958 

0.952 

0.964 

6 

45 

0.476 

0.529 

0.434 

0.623 

7 

-45 

0.688 

0.717 

0.623 

0.811 

8 

90 

0.993 

0.994 

0.988 

1.000 

9 

90 

0.993 

0.994 

1.000 

0.988 

10 

-45 

0.675 

0.717 

0.811 

0.623 

11 

45 

0.463 

0.529 

0.623 

0.434 

12 

0 

0.953 

0.958 

0.964 

0.952 

13 

90 

0.940 

0.946 

0.952 

0.939 

14 

-45 

-0.17 

-0.037 

0.057 

-0.131 

15 

45 

-0.39 

-0.226 

-0.132 

-0.320 

16 

0 

0.900 

0.909 

0.915 

0.903 
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Figure 4.8: Graph of M12 for Cyclic in Plane Bending Load 

From the figures we observe that the plies (0/45/-45/90/0) from the bottom and (0/90/- 
45/45/0) top fail in case of application of bending loads in the longitudinal direction, which is 
expected while the transverse direction has the plies (90/-45/45/0) from the top failing. In case of 
shear the plies (-45/45) fail from the ends followed by the (-45/45) plies in the center. Also we 
observe that in case of all the predictions for slice 1 and slice 125 the ICAN/JAVA with no sub- 
slicing predicts the average of the two values, which is consistent throughout. 

4.6.3 MFIM Results 

The MFIM is a new addition to ICAN/JAVA and a durability analysis was carried out to 
observe the effects of MFIM on the results. Figures 4.9-4.11 show the plots for the various 
factors of safety. 
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Figure 4.9: Graph of Ml 1 for MFIM 



Figure 4.10: Graph of M22 for MFIM 
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1 



1(0 DEG) 2(90 DEG) 3 (90 DEG) 4 (0 DEG) 

PLY-LAYUP 

Figure 4.1 1: Graph of M12 for MFIM 

In Figure 5.16, the stresses in the 11 directions vary in their values by only 5%. The 
difference is reduced to about 0.5% in Figure 5.17 for the case of stresses in 22 directions. While 
the shear stresses show a difference of 16% in the center of the ply layup. This shows that MFIM 
indeed has some effect on the results for the same models and loads. An attempt was made to try 
to bring the two variations as close by using different coefficients for MFIM modeling. We can 
see that the 0 degree plies fail first in shear and transverse directions while in case of longitudinal 
direction the 0 degree ply in the bottom fails first. 

4.7 Conclusion 

From the above comparisons we can observe that there are some differences in the results of the 
software ICAN and ICAN/JAVA. This may be due to the addition of new modules like sub- 
slicing, variable loading etc in ICAN/JAVA. The other cause can be some bugs that need to be 
fixed. We also observe that the MFIM has an effect on the results as the material properties are 
changed to a considerable extent. This is important because then we can take into account the 
change in material properties with application of loads over time to get more accurate results. 
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4.8 Future Work 


Hence, a further detailed look into the new ICAN/JAVA for its corroboration with the ICAN 
results is required to fix any bugs that might be present in the new software. Also a need to 
further study in detail the new enhancements of ICAN/JAVA for better results is necessary to 
realize the full potential of the software. 
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Chapter 5 

Prediction of First Ply Failure and Fracture in Composite 
Materials of Different Size and Geometry 


5.1 Objective 

The objective of this chapter is to simulate the biaxial failure of composites of different size 
and geometry and compare them with experimental results. The research consists of generation 
of first-ply failure envelopes for combined loading of these laminates on the basis of first-ply 
failure and laminate fracture. 

5.2 Background 

For the last three decades or so, efforts have been made to predict the laminate fracture under 
uniaxial and/or combined loading. Due to the existence of a large number of failure criteria and 
analysis methods no exact solution has been found which can correctly predict the composite 
laminate fracture. But a comparison of the various failure theories with a benchmarked 
experimental data can be useful and instructive for structural designers. 

To that end, Hinton and Soden (27) organized an exercise to compare the current theories of 
composite failure with their experimental results. Among the many theories Gotsis et al. (28) 
used the micromechanics-based theory and progressive fracture. The results were comparable 
with the experimental data but not very accurate for the biaxial failure envelopes. Others like Sun 
et al. (29) used the linear laminate analysis, Rotem used the Rotem failure criteria (30) and 
Worlfe et al. (31) used the strain-energy based criteria to compare with the experimental data. 
Some of the theories could accurately predict the stress/strain curves while others could predict 
the biaxial failure but no one theory could predict the laminate fracture for all the benchmarks 
and hence only a comparison of many theories occurred. An attempt has been made to predict 
the biaxial failure envelope for composites using the CODSTRAN and ICAN in this research and 
to compare with the data of Hinton et al. (27) 
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5.3 Theory 


CODSTRAN is used to simulate multi-layered fiber composite structures under any loading 
and hygrothennal conditions. It pennits a simulation of local behavior as well as global structural 
behavior as shown in Figure 5.1 



Figure 5.1: Simulation Cycle of Progressive Fracture in Composite 
Laminates and Structures via CODSTRAN 

The constituent properties are described as functions of environmental and mechanical 
loading conditions as shown in the lower part of the figure. By knowing the developed ply 
stresses, and ply strength we can predict the damage initiation, growth, accumulation, and 
propagation, which results in the constituent material properties being updated at every load 
increment. For example, if a ply’s transverse strength exceeds its allowable strength then the ply 
is assumed to fail and the matrix modulus is updated as negligible. But if a ply longitudinal stress 
exceeds the allowable strength, then the fiber/matrix longitudinal moduli are both updated as 
negligible. This means that the ply is assumed to carry no load and the stresses are then 
redistributed to the surrounding plies. Once the current constituent properties are determined the 
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repeated application of micromechanics, macromechanics, and laminate theory are used to 
assemble the global structural stiffness matrix, which is fed to the finite element analysis. 

The nodal stress resultants are from the finite element analysis and are used to decompose the 
changes in the global response (laminate stress and strains) on the local material stress/resistance. 
The load is increased only if there is no further damage due to changes in ply level stresses. 
Otherwise, only the material properties are updated at every increment until equilibrium is 
reached between the applied load and the local response. Until a specified convergence is 
reached the structural equilibrium is maintained as shown in Figure 5.2. 



Figure 5.2: CODSTRAN Load Increment 

The final result in terms of load as a function of global displacement is shown in Figure 5.3. 
Agarwal et al. (32) have described the damage progression analysis procedure in detail. 
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Figure 5.3: Overall CODSTRAN Simulations 


5.4 Model Description 

Two different specimens were tested under biaxial stress state to generate the failure 
envelopes. The material used was E-glass/MY750 [±55] s (Appendix B) laminate. A square one- 
inch by one-inch coupon was tested under transverse and shear loading as shown in Figure 5.4. 
The coupon model was similar to the one used by Gotsis et al. (28) to represent the stress state 
prior to the disclosure of the cylindrical specimen and it’s boundary conditions, by Hinton and 
Soden (27). 
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The results were further compared with a cylindrical laminate subjected to a uniform internal 
pressure. The length of the specimen is 8.4 inches and the tube is clamped at both ends. The 
inner diameter of the tube is 39.4 inches with a thickness of 0.4 inches as shown in Figure 5.5. 
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Figure 5.5: E-glass/MY750 [±55] s Cylindrical Tube 


End reinforcement with variable thickness made of E-glass/MY750 [90] is applied to the 
tube. The reinforcement is provided to simulate real testing conditions where the specimen is 
held between the ends. The main focus is at the center of the tube for a length of 23.6 inches 
where the E-glass/MY750 [±55] s laminate is present with no reinforcement. There laminate 
layups used the base of [±55] s as the basic core and the reinforcements were added as 90° lamina 
each of thic kn ess 0.01578 inch. To simulate the variable reinforcement twelve ply layup were 
defined and were applied to different sections of the cylinder as shown in the Figure 5.6. Hence, 
the thickness of the cylinder varied from 0.039 inches at the center to 0.25 at the ends as shown 
in Figure 5.7. 
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Figure 5.6: Twelve Ply Schedules to Simulate Reinforcement for Cylinder 
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Figure 5.7: Variable Thickness of the Cylinder Due to Different Ply-Schedules 
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The cylindrical tube was further tested with the calibrated properties SGE1/MHD1 and the 
results compared with the non-calibrated material properties. These material properties are 
shown in Appendix B. 

5.5 Model Assumptions 

The coupon was fixed in the z direction along the boundary. The first node at the origin was 
fixed in all directions while the x-axis and y-axis were fixed in y and x directions respectively. 
The thickness of the coupon was 0.039 inches, which is the same as the thickness of the center 
part of the cylinder. Axial compression and tension were applied to the coupon along the x and y 
directions to simulate the stress state. A force of 1 kip was used as the initial load increment. 
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Figure 5.8: Boundary Conditions for the Coupon 
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To simulate the boundary conditions of clamped specimen for the cylinder one end of it was 
fixed in all directions while the other end was fixed in y and z directions only as shown in Figure 
5.8. This allowed for deformation of the specimen in the axial direction. 



Figure 5.9: E-glass/MY750 Tube Geometry 


The axial force and internal pressure were applied as loads to simulate the experimental 
conditions as shown in Figure 5.9 and the results plotted. 
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Figure 5.10: Direction of Force and Pressure Applied to the Cylinder Specimen 


The stress ratio was varied from -1 to 1 and the following procedure used to calculate the 
required internal pressure. 


°y = 


pr 

t 

SC 

2TL rt 


StressRatio = 




(5.1) 


where, 

p = int ernal _ pressure 
r = radius _of _ cyclinder 
t = thickness _of _ cyclinder 


The force applied was 1 kip and the total force on the top of the cylinder was 32 kips. We 
calculated the internal pressure for SR=2/1 as follows. 
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2'XA 

2Y\rtr 


fir 2 


32 

n(i.97) 2 


2.62 psi / inch 2 


(5.2) 


where, 

p = pressure 

y F x = summation _of _ forces 
r = radius of _cyclinder = 1. 97 inch 
t = thickness _ of _cyclinder = 0.039 inch 

5.6 Material Calibration 

As CODSTRAN is micromechanics based code, a calibration of the material properties so as 
to match the predicted composite unidirectional properties were undertaken. The following graph 
shown in Figure 5.10 was used to calibrate the properties. 



Figure 5.11: Biaxial Envelope of 0° GRP Lamina Under 
Combined c x and c y Stresses for E-glass/MY750 Epoxy. 
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ICAN was used to carry out the material calibration as it has the ability to calculate the 
laminate properties from the fiber and matrix properties. The model used for the calibration was 
a four ply 0° unidirectional laminate. The existing material properties of SGES/MHDY 
(Appendix B) were used as the basis of our calibration. In-plane loads were applied to check the 
response to the variation of the material properties. 

Table 5.1 shows the comparison of the stresses obtained from simulated and corrected 
material properties to that of experimental data. 


Table 5.1: Experimental and simulated data for material calibration 


PMEMB 

Stress (MPA) Experimental 

Stress (MPA) Simulated 

In x-direction 

In y-direction 

Sigma x 

Sigma y 

Sigma x 

Sigma y 

+200 

0 

+1271.86 

0.000 

+1290.545 

0.000 

-200 

0 

-830.508 

0.000 

-773.619 

0.000 

0 

+10 

0.000 

+42.373 

0.000 

+66.647 

0 

-10 

0.000 

-140.678 

0.000 

-99.975 


From Figure 5.11 we observe the difference in the biaxial failure envelope for the lamina to 
be less than 7% in the sigma x direction and less than 40% in the sigma y direction with respect 
to the experimental data. 
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—♦—Experimental —■— Simulated 



Figure 5.12: Comparison of Biaxial Envelope of 0° GRP Lamina Under Combined Ox and Cy 
Stresses for Experimental and Simulated Material Properties. 

The error in the y direction was higher because a slight change in properties could increase 
the error of the x direction and this was the closest we could get to match the envelope overall. 
Tables 5.2 and 5.3 show the new material properties and Chapter 5.10 has the original properties. 
Comparisons of the fabric and matrix properties are shown in Table 5.4 and 5.5. 
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Table 5.2: SGE1 -Fiber Properties 

SGE1 Silenka E-glass 1200 tex Fiber (Soden). 

$ 

$ 

$ 


Number of fibers per end 

Nf 

3000 

number 

Filament equivalent diameter 

df 

0.300E-03 

inches 

Weight density 

Rhof 

0.640E-01 

lb/in**3 

Normal moduli ( 1 1) 

Efl 1 

1.073E+07 

psi 

Normal moduli (22) 

Ef22 

1.073E+07 

psi 

Poisson’s ratio (12) 

Nufl2 

0.200E+00 

non-dim 

Poisson’s ratio (23) 

Nuf23 

0.200E+00 

non-dim 

Shear moduli (12) 

Gfl2 

4.467E+06 

psi 

Shear moduli (23) 

Gf23 

4.467E+06 

psi 

Thermal expansion coef. (11) 

Alfafl 1 

0.272E-05 

in/in/F 

Thermal expansion coef. (22) 

Alfaf22 

0.272E-05 

in/in/F 

Heat conductivity (11) 

Kfl 1 

5.208E-02 

BTU-in/hr/in 2 /F 

Heat conductivity (22) 

Kf22 

5.208E-02 

BTU-in/hr/in 2 /F 

Heat capacity 

Cf 

0.170E+00 

BTU/lb/F 

Fiber tensile strength 

SfT 

3.073E+05 

psi 

Fiber compressive strength 

SfC 

2.103E+05 

psi 



Table 5.3: 

MHDl -Matrix Properties 

MHDl My750/HY 9 1 7/D Y 063 Matrix (Soden) 
$ 



$ 

Weight density 

Rhom 

0.470E-01 

lb/in**3 

Normal modulus 

Em 

0.486E+06 

psi 

Poisson"s ratio 

Num 

0.350E+00 

non-dim 

Thermal expansion coef. 

Alfa m 

0.322E-04 

in/in/F 

Matrix heat conductivity 

Km 

8.681E-03 

BTU-in/hr/in 2 /F 

Heat capacity 

Cm 

0.250E+00 

BTU/lb/F 

Matrix tensile strength 

SmT 

0.737E+04 

psi 

Matrix compressive strength 

SmC 

0.245E+05 

psi 

Matrix shear strength 

SmS 

0.122E+05 

psi 

Allowable tensile strain 

eps mT 

0.500E-01 

in/in 

Allowable compr. strain 

eps mC 

0.500E-01 

in/in 

Allowable shear strain 

eps mS 

0.500E-01 

in/in 

Allowable torsional strain 

eps mTOR 

0.500E-01 

in/in 

Void heat conductivity 

kv 

0.225E+00 

BTU-in/hr/in 2 /F 

Glass transition temperature 

Tgdr 

0.420E+03 

F 
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Table 5.4: Comparison of Fiber Properties with Respect to Reference Material 


Properties 

Reference Material 

Calibrated Material 

Normal Moduli Ell psi 

1.073E+07 

1.073E+07 

Normal Moduli E22 psi 

1.073E+07 

1.073E+07 

Poisson’s Ratio NU12 

0.200E+00 

0.200E+00 

Poisson’s Ratio NU23 

0.200E+00 

0.200E+00 

Shear Moduli G12 psi 

4.467E+06 

4.467E+06 

Shear Moduli G23 psi 

4.467E+06 

4.467E+06 

Tensile Strength ST psi 

3.118E+05 

3.073E+05 

Compressive Strength SC psi 

2.103E+05 

2.103E+05 


Table 5.5: Comparison of Matrix Properties with Respect to Reference Material 


Properties 

Reference material 

Calibrated material 

Normal Modulus E psi 

0.486E+06 

0.486E+06 

Poisson’s Ratio NU 

0.3500E+00 

0.3500E+00 

Tensile Strength ST psi 

0.149E+05 

0.737E+04 

Compressive Strength SC psi 

0.174E+05 

0.245E+05 

shear strength SS psi 

0.100E+05 

0.122E+05 


5.7 Results and Discussion 


The simulations were done in three parts. Fist the various uniaxial and biaxial loads were 
applied to the coupon and the results for the first-ply failure were converted to the stresses in x 
and y directions and tabulated in Table 5.6. 


Table 5.6: Coupon simulation Data 


Applied Initial Force 
Increment (kip) 

Resultant Failure 
Force (kip) 

Failure 

(M 

■ Stress 
pa) 

In X-Direction 

In Y-Direction 

In X-Direction 

In Y-direction 

Sigma X 

Sigma Y 

+1 

0 

+0.4429 

0.0000 

+078.30 

0.00 

-1 

0 

-0.3907 

0.0000 

-69.07 

0.00 

0 

+1 

0.0000 

+1.4770 

0.00 

+261.10 

0 

-1 

0.0000 

-3.1630 

0.00 

-559.18 

-1 

-1 

-1.9560 

-1.9560 

-345.80 

-345.80 

+1 

+1 

+1.6540 

+1.6540 

+292.40 

+292.40 

+1 

-1 

+2.2310 

-2.2310 

+394.40 

-394.40 

-1 

+1 

-0.4387 

+0.4387 

-077.55 

+077.55 
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The following fonnula was used to convert the force into stress for the coupon. 


y F x 

a = 6.8941 5Mpa 

where, 

L = length _ of _ coupon = 1.0 inch 
t = thickness _of _coupon = 0.039 inch 


(5.3) 


Secondly the non-calibrated properties of SGES/MHDY laminate were used to simulate the 
biaxial/uniaxial loading and tabulated in Table 5.7. 


Table 5.7: Cylinder Simulation Data: SGES/MHDY Laminate 


Initial Force 
Increment in 
X-Direction 
(kip) 

Applied 

Internal 

Pressure 

(psi) 

Resultant 
Failure Force 
in X-Direction 
(kip) 

Resultant 

Failure 

Pressure 

(psi) 

Sigma X 
(Mpa) 

Sigma Y 
(Mpa) 

+1 

0.0000 

+5.842 

0.0000 

+83.44 

0.00 

-1 

0.0000 

-6.615 

0.0000 

-94.98 

0.00 

+1 

-1.3123 

+3.655 

-0.1499 

+52.20 

-52.20 

-1 

+1.3123 

-4.513 

+0.1851 

-64.46 

+64.46 

+1 

+1.3123 

+7.179 

+0.2944 

+102.54 

+102.53 

-1 

-1.3123 

-9.725 

-0.3988 

-138.91 

-138.89 


+1.3123 

0.000 

+0.5024 

0.00 

+174.97 


-1.3123 

0.000 

-0.4061 

0.00 

-141.43 

+1 

+2.6246 

+10.930 

+0.8961 

+156.12 

+312.08 

+1 

+2.2730 

+11.830 

+0.8400 

+168.97 

+292.55 

+1 

+4.5931 

+7.280 

+1.0450 

+103.98 

+363.94 


Thirdly the calibrated properties were subjected to the same loads and the results tabulated in 
Table 5.8 
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Table 5.8: Cylinder Simulation Data: SGE1/MHD1 Laminate 


Initial Force 
Increment in 
X-Direction 
(kip) 

Applied 

Internal 

Pressure 

(psi) 

Resultant 
Failure Force 
in X-Direction 
(kip) 

Resultant 

Failure 

Pressure 

(psi) 

Sigma X 
(Mpa) 

Sigma Y 
(Mpa) 

+1 

0.0000 

+4.565 

0.0000 

+65.21 

0.00 

-1 

0.0000 

-8.379 

0.0000 

-119.68 

0.00 


+1.3123 

0.000 

+0.5450 

0.00 

+189.80 


-1.3123 

0.000 

-0.3024 

0.00 

-105.32 

+1 

+1.3123 

+5.190 

+0.2129 

+74.13 

+74.14 

+1 

-1.3123 

+2.513 

-0.1030 

+35.89 

-35.87 

-1 

+1.3123 

-5.469 

+0.2243 

-78.11 

+78.11 

-1 

-1.3123 

-13.720 

-0.5626 

-195.97 

-195.94 

-1 

-3.2808 

-7.744 

-0.7940 

-110.61 

-276.53 

+1 

+3.2808 

+7.929 

+0.8129 

+113.25 

+283.11 

+1 

+4.2808 

+6.833 

+0.9141 

+97.60 

+318.35 

+1 

+8.2808 

+2.288 

+0.5920 

+32.68 

+206.17 

+1 

+0.2808 

+4.727 

+0.0415 

+67.52 

+14.440 

+1 

+4.2809 

+6.844 

+0.9155 

+97.75 

+318.84 

+1 

+4.2812 

+6.833 

+0.9142 

+97.60 

+318.39 

+1 

+4.3000 

+5.783 

+0.7771 

+82.60 

+270.64 

+1 

+4.2830 

+6.703 

+0.8971 

+95.74 

+312.43 

+1 

+4.2815 

+6.740 

+0.9018 

+96.27 

+314.07 

+1 

+4.2811 

+6.182 

+0.8270 

+88.30 

+288.02 

+1 

+4.2812 

+6.833 

+0.9142 

+97.60 

+318.39 

+1 

+4.2800 

+6.029 

+0.8063 

+86.11 

+280.81 

+1 

+4.2810 

+6.844 

+0.9156 

+97.75 

+318.88 


The following fonnula was used to convert the force into stress for the cylinder. 


°x 


Yf, 

— — — * 6.89475Mpa 
m rt 


(5. 


where, 

r = radius _ of _cyclinder = 1.97 inch 
t = thickness _of _ cyclinder = 0.039 inch 


The following formula was used to convert the pressure into Mega Pascal. 
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— * 6.89415Mpa 
t 


(5.5) 


cr y = 


Figure 5.12 and 5.13 show the comparison of the experimental and theoretical results for the 
biaxial failure envelopes of E-glass/MY750 [±55] s laminate, coupon and the SGE1/MHD1 
[±55] s laminate under combined a x - a y biaxial loading. 


—♦—Coupon (SGES/MHDY) —■—Cylinder (SGES/MHDY) Cylinder (SGE1/MHD1)] 



Figure 5.13: Biaxial Failure Stresses for [±55] s E-glass/MY750 
Laminate Using CODSTRAN 
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Figure 5.14: Biaxial Failure Stresses for [±55] sE-glass/MY750 
Laminate Obtained Experimentally 


The theoretical final strength computed by CODSTRAN under predicted the experimental 
leakage strength between SR=1:1 and SR=2:1. In the compression-compression quadrant 
CODSTRAN predicted that the composite had no damage tolerance or the ultimate fracture to 
occur immediately after first ply failure due to longitudinal compressive stresses. These 
predictions were approximately three times lower than the test results at a stress ratio of 
approximately -2:1. This may be attributed to (a) the high fiber volume fraction of the test 
specimens used in this quadrant [16], (b) the use of thick tubes in the experiments while the 
theory uses thin flat plates, (c) simulations neglected the residual matrix stiffness of the failed 
plies that may have under predicted the ultimate failure load under compressive stresses. The 
final strength for the coupon as predicted by CODSTRAN is higher than the cylindrical 
specimen because the coupon is not subjected to internal pressure but a combination of o x - a y 
stresses. The additional boundary conditions used alongside the edges of the coupon, were 
responsible for the stiffening of the coupon resulting in higher values for ultimate failure load. 
There were no experiments for the coupon and only the cylinder results were compared with the 
given experimental data. 
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5.8 Conclusion 


A comparison was made between the theoretical predictions and progressive fracture 
simulations using CODSTRAN and it was found that only a limited number of segments of the 
biaxial failure envelope for [±55] s glass/epoxy laminate correlated to the test results. The final 
strengths were under-predicted in a vast majority of segments. This can be attributed to the high 
fiber volume fraction of the test specimens, the use of thick tubes in the experiments while the 
simulations use a thin plate and the simulations neglected the residual matrix stiffness of the 
failed plies that may have predicted the ultimate failure loads under compressive stresses. 

The ability of CODSTRAN to predict first ply failure envelopes is fairly accurate provided 
the structural condition of the specimen is known. 

5.9 Future Work 

A need to simulate the exact structural conditions at the time of experimentation to compare 
with the experimental results is required. The different models like thin and thick tubes need to 
be taken into account in the simulation to get more accurate comparison. Residual matrix 
stiffness after transverse tensile failure needs to be quantified and accounted for in the damage 
simulations. 
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5.10 Material Properties 


Fiber Properties: 


AS-4 

Number of fibers per end = 10000 
Fiber diameter = 0.00762 mm (0.300E-3 in) 

Fiber Density = 4.04E-7 Kg/m 3 (0.063 lb/in 3 ) 

Longitudinal normal modulus = 227 GPa (32.90E+6 psi) 

Transverse normal modulus = 13.7 GPa (1.99E+6 psi) 

Poisson's ratio (v 12 ) = 0.20 
Poisson's ratio (v 23 ) = 0.25 
Shear modulus (Gi 2 ) = 13.8 GPa (2.00E+6 psi) 

Shear modulus (G 23 ) = 6.90 GPa (1.00E+6 psi) 

Longitudinal thermal expansion coefficient = -1.0E-6/°C (-0.55E-6 /°F) 
Transverse thermal expansion coefficient = 1.0E-5/°C (0.56E-5 /°F) 
Longitudinal heat conductivity = 43.4 J-m/hr/m 2 /°C (580 BTU-in/hr/in 2 /°F) 
Transverse heat conductivity = 4.34 J-m/hr/m 2 /°C (58 BTU-in/hr/in 2 /°F) 
Heat capacity = 0.712 KJ/Kg/°C (0.17 BTU/lb/°F) 

Tensile strength = 3.723 GPa (540 ksi) 

Compressive strength = 3.351 GPa (486 ksi) 


SGES Silenka E-glass 1200 tex Fiber (Soden). 


Number of fibers per end 

Nf 

3000 

number 

Filament equivalent diameter 

df 

0.300E-03 

inches 

Weight density 

Rhof 

0.640E-01 

lb/in**3 

Normal moduli ( 1 1) 

Efl 1 

1.073E+07 

psi 

Normal moduli (22) 

Ef22 

1.073E+07 

psi 

Poisson"s ratio (12) 

Nufl2 

0.200E+00 

non-dim 

Poisson"s ratio (23) 

Nuf23 

0.200E+00 

non-dim 

Shear moduli (12) 

Gfl2 

4.467E+06 

psi 

Shear moduli (23) 

Gf23 

4.467E+06 

psi 

Thermal expansion coef. (11) 

Alfafl 1 

0.272E-05 

in/in/F 

Thermal expansion coef. (22) 

Alfaf22 

0.272E-05 

in/in/F 

Heat conductivity (11) 

Kfl 1 

5.208E-02 

BTU-in/hr/in**2/F 

Heat conductivity (22) 

Kf22 

5.208E-02 

BTU-in/hr/in**2/F 

Heat capacity 

Cf 

0.170E+00 

BTU/lb/F 

Fiber tensile strength 

SfT 

3.118E+05 

psi 

Fiber compressive strength 

SfC 

2.103E+05 

psi 


T300 Graphite Fiber. 

Number of fibers per end 

Nf 

3000 

number 

Filament equivalent diameter 

df 

0.300E-03 

inches 

Weight density 

Rhof 

0.640E-01 

lb/in**3 

Normal moduli ( 1 1) 

Efl 1 

0.320E+08 

psi 

Normal moduli (22) 

Ef22 

0.200E+07 

psi 

Poisson"s ratio (12) 

Nufl2 

0.200E+00 

non-dim 

Poisson"s ratio (23) 

Nuf23 

0.250E+00 

non-dim 

Shear moduli (12) 

Gfl2 

0.130E+07 

psi 
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Shear moduli (23) 

Gf23 

0.700E+06 

psi 

Thermal expansion coef. (11) 

Alfafl 1 

-.550E-06 

in/in/F 

Thermal expansion coef. (22) 

Alfaf22 

0.560E-05 

in/in/F 

Heat conductivity (11) 

Kfl 1 

0.403E+01 

BTU-in/hr/in* *2/F 

Heat conductivity (22) 

KT22 

0.403E+00 

BTU-in/hr/in**2/F 

Heat capacity 

Cf 

0.170E+00 

BTU/lb/F 

Fiber tensile strength 

SfT 

0.350E+06 

psi 

Fiber compressive strength 

SfC 

0.300E+06 

psi 


AS-- Graphite Fiber. 


Number of fibers per end 

Nf 

10000 

number 

Filament equivalent diameter 

df 

0.300E-03 

inches 

Weight density 

Rhof 

0.630E-01 

lb/in* *3 

Normal moduli ( 1 1) 

Efl 1 

0.310E+08 

psi 

Normal moduli (22) 

Ef22 

0.200E+07 

psi 

Poisson"s ratio (12) 

Nufl2 

0.200E+00 

non-dim 

Poisson"s ratio (23) 

Nuf23 

0.250E+00 

non-dim 

Shear moduli (12) 

Gfl2 

0.200E+07 

psi 

Shear moduli (23) 

Gf23 

0.100E+07 

psi 

Thermal expansion coef. (11) 

Alfafl 1 

-.550E-06 

n/in/F 

Thermal expansion coef. (22) 

Alfaf22 

0.560E-05 

in/in/F 

Heat conductivity (11) 

Kfl 1 

0.403E+01 

BTU-in/hr/in**2/F 

Heat conductivity (22) 

Kf22 

0.403E+00 

BTU-in/hr/in**2/F 

Heat capacity 

Cf 

0.170E+00 

BTU/lb/F 

Fiber tensile strength 

SfT 

0.400E+06 

psi 

Fiber compressive strength 

SfC 

0.400E+06 

psi 


SGLA S-Glass Fiber. 


Number of fibers per end 

Nf 

204 

number 

Filament equivalent diameter 

df 

0.360E-03 

inches 

Weight density 

Rhof 

0.900E-01 

lb/in* *3 

Normal moduli ( 1 1) 

Efll 

0.124E+08 

psi 

Normal moduli (22) 

Ef22 

0.124E+08 

psi 

Poisson"s ratio (12) 

Nufl2 

0.200E+00 

non-dim 

Poisson"s ratio (23) 

Nuf23 

0.200E+00 

non-dim 

Shear moduli (12) 

Gfl2 

0.517E+07 

psi 

Shear moduli (23) 

Gf23 

0.517E+07 

psi 

Thermal expansion coef. (11) 

Alfafl 1 

0.280E-05 

in/in/F 

Thermal expansion coef. (22) 

Alfaf22 

0.280E-05 

in/in/F 

Heat conductivity (11) 

Kfl 1 

5.208E-02 

BTU-in/hr/in**2/F 

Heat conductivity (22) 

Kf22 

5.208E-02 

BTU-in/hr/in**2/F 

Heat capacity 

Cf 

0.170E+00 

BTU/lb/F 

Fiber tensile strength 

SfT 

0.360E+06 

psi 

Fiber compressive strength 

SfC 

0.300E+06 

psi 
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Matrix Properties: 


5250 HMHS 

Matrix density = 3.40E-7 Kg/m 3 (0.0457 lb/in 3 ) 

Normal modulus = 4.62 GPa (671 ksi) 

Poisson's ratio = 0.705 

Coefficient of thermal expansion = 0.518E-4/°C (0.288E-4 /°F) 

Heat conductivity = 0.649E-3 J-m/hr/m 2 /°C (0.862E-8 BTU-in/hr/in 2 /°F) 
Heat capacity = 1.047 KJ/Kg/°C (0.25 BTU/lb/°F) 

Tensile strength = 75.1 MPa (10.9 ksi) 

Compressive strength = 283 MPa (41.0 ksi) 

Shear strength =138 MPa (20.0 ksi) 

Allowable tensile strain = 0.02 
Allowable compressive strain = 0.05 
Allowable shear strain = 0.04 
Allowable torsional strain = 0.04 

Void conductivity = 16.8 J-m/hr/m 2 /°C (0.225 BTU-in/hr/in 2 /°F) 

Glass transition temperature = 216°C (420°F) 

MHDY My750/HY917/DY063 Matrix (Soden) 


Weight density 

Rhom 

0.470E-01 

lb/in**3 

Normal modulus 

Em 

0.486E+06 

psi 

Poisson"s ratio 

Num 

0.350E+00 

non-dim 

Thermal expansion coef. 

Alfa m 

0.322E-04 

in/in/F 

Matrix heat conductivity 

Km 

8.681E-03 

BTU -in/hr/in* *2/F 

Heat capacity 

Cm 

0.250E+00 

BTU/lb/F 

Matrix tensile strength 

SmT 

0.149E+05 

psi 

Matrix compressive strength 

SmC 

0.174E+05 

psi 

Matrix shear strength 

SmS 

0.100E+05 

psi 

Allowable tensile strain 

eps mT 

0.500E-01 

in/in 

Allowable compr. strain 

eps mC 

0.500E-01 

in/in 

Allowable shear strain 

eps mS 

0.500E-01 

in/in 

Allowable torsional strain 

eps mTOR 

0.500E-01 

in/in 

Void heat conductivity 

kv 

0.225E+00 

BTU -in/hr/in* *2/F 

Glass transition temperature 

Tgdr 

0.420E+03 

F 

IMLS Intermediate Modulus Low Strength Matrix. 



Weight density 

Rhom 

0.460E-01 

lb/in**3 

Normal modulus 

Em 

0.500E+06 

psi 

Poisson"s ratio 

Num 

0.410E+00 

non-dim 

Thermal expansion coef. 

Alfa m 

0.570E-04 

in/in/F 

Matrix heat conductivity 

Km 

8.681E-03 

BTU-in/hr/in* *2/F 

Heat capacity 

Cm 

0.250E+00 

BTU/lb/F 

Matrix tensile strength 

SmT 

0.700E+04 

psi 

Matrix compressive strength 

SmC 

0.210E+05 

psi 

Matrix shear strength 

SmS 

0.700E+04 

psi 

Allowable tensile strain 

eps mT 

0.140E-01 

in/in 

Allowable compr. Strain 

eps mC 

0.420E-01 

in/in 

Allowable shear strain 

eps mS 

0.320E-01 

in/in 

Allowable torsional strain 

eps mTOR 

0.320E-01 

in/in 

Void heat conductivity 

kv 

0.225E+00 

BTU -in/hr/in* *2/F 

Glass transition temperature 

Tgdr 

0.420E+03 

F 
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IMHS Intermediate Modulus High Strength Matrix. 


Weight density 

Rhom 

0.440E-01 

lb/in**3 

Normal modulus 

Em 

0.500E+06 

psi 

Poisson"s ratio 

Num 

0.350E+00 

non-dim 

Thermal expansion coef. 

Alfa m 

0.360E-04 

in/in/F 

Matrix heat conductivity 

Km 

8.681E-03 

BTU-in/hr/in* *2/F 

Heat capacity 

Cm 

0.250E+00 

BTU/lb/F 

Matrix tensile strength 

SmT 

0.150E+05 

psi 

Matrix compressive strength 

SmC 

0.350E+05 

psi 

Matrix shear strength 

SmS 

0.130E+05 

psi 

Allowable tensile strain 

eps mT 

0.200E-01 

in/in 

Allowable compr. strain 

eps mC 

0.500E-01 

in/in 

Allowable shear strain 

eps mS 

0.350E-01 

in/in 

Allowable torsional strain 

eps mTOR 

0.350E-01 

in/in 

Void heat conductivity 

kv 

0.225E+00 

BTU -in/hr/in* *2/F 

Glass transition temperature 

Tgdr 

0.420E+03 

F 


HMHS High Modulus High Strength Matrix. 

$ 

$ 

$ 


Weight density 

Rhom 

0.450E-01 

lb/in**3 

Normal modulus 

Em 

0.750E+06 

psi 

Poisson"s ratio 

Num 

0.350E+00 

non-dim 

Thermal expansion coef. 

Alfa m 

0.400E-04 

in/in/F 

Matrix heat conductivity 

Km 

8.681E-03 

BTU -in/hr/in* *2/F 

Heat capacity 

Cm 

0.250E+00 

BTU/lb/F 

Matrix tensile strength 

SmT 

0.200E+05 

psi 

Matrix compressive strength 

SmC 

0.500E+05 

psi 

Matrix shear strength 

SmS 

0.150E+05 

psi 

Allowable tensile strain 

eps mT 

0.200E-01 

in/in 

Allowable compr. strain 

eps mC 

0.500E-01 

in/in 

Allowable shear strain 

eps mS 

0.400E-01 

in/in 

Allowable torsional strain 

eps mTOR 

0.400E-01 

in/in 

Void heat conductivity 

kv 

0.225E+00 

BTU -in/hr/in* *2/F 

Glass transition temperature 

Tgdr 

0.420E+03 

F 


NASA/CR— 2003-212351 


97 



5.11 References 


1. Chamis, C.C: Computer Code for the Analysis of Multilayered Fiber Composites-User’s 
Manual. NASA TN-D-70 13, 1971. 

2. Chamis, C.C., and Sinclair, J.H.: Computer Code for Intraply Hybrid Composite Design. 
NASA Tm-82593, 1981. 

3. Chamis, C.C., and Sinclair, J.H.: INHYD: Computer Code for Intraply Hybrid Composite 
Design. NASA TP-2239, 1983. 

4. Murthy, P.L.N., and Chamis, C.C.: Integrated Composite Analyzer ICAN: User’s And 
Programmer’s Manual. NASA Tp-2512, 1986. 

5. Chamis, C.C., and Smith, G.T.: Composite Durability Structural Analysis. NASA TM- 
79070, 1978. 

6. Nakazawa, S., Dias, J.B., and Spiegel, M.S.: MHOST User’s Manual. NASA-Marc Analysis 
Research Corporation, 1987. 

7. CSTEM User’s Manual, NASA NAS3-26385, 2000 

8. Chamis, C.C. Tntegrated Composite Analyzer (ICAN/Java)-User’s Manual. NASA, 2001. 

9. Rubenstien, Robert, Chamis, C. C.: STAEBL/General Composites With Hydrothermal 
Effects: NASA Technical Memorandum, Dec, 1987, lOp 

10. Vanderplaats, G.N.: CONMIN, A Fortran Program For Constrained Function Minimization. 
User's Manual. NASA Technical Memorandum, Tm— X— 62, 282, 1973. 

11. Alpha Star Corporation: Genoa/Progressive Failure Analysis Software, Long Beach, Ca 

12. Anderberg, Y: User’s Manual For TCD 3.0 With Super- TEMPCALC, Fire Safety Design, 
Lund, Sweden, 1990 

13. Liu, T.C.H: Three-Dimensional Modeling Of Steel/Concrete Composite Connection 
Behavior In Fire, Journal Of Construction Steel Research, V46, Apr-Jun 1998 

14. Bailey, C.G, Burgess, I.W., Plank, R.J.: Computer Simulation Of Full-Scale Structural Fire 
Test, Structural Engineer, V74, Mar 1996 

15. Lamont, S., Usmani, A.S., Drysdale, D.D.: Heat Transfer Analysis Of The Composite Slab In 
The Cardington Frame Fire Tests, Fire Safety Journal, V3, Nov 2001 

16. Looyeh, M.R.E., Bettess, P, Gibson, A.G.: One-Dimensional Finite Element Simulation For 
The Fire-Performance Of GRP Panels For Offshore Structures, Int. Journal Of Numerical 
Methods For Heat And Fluid Flow, V7, 1997, Pg 609-625 

17. Asaro, Robert J, Dao, Ming: Fire Degradation Of Fiber Composites, Marine Technology, 
V34,Jul 1997, Pg 197-210 

18. Dembsey, N. A., Jacoby, D.J.: Evaluation Of Common Ignition Models For Use With Marine 
Cored Composites, Fire And Materials, V24, 2000 

19. Kreith, F: CRC Handbook Of Thermal Engineering, CRC Press, Dec 1999 

20. Lattimer, B.Y, Fire Resistance Testing Of Flat Horizontal Sandwich Panel With Load, 
Hughes Associates, May 2002 

21. Yang, F., Li, J. C. M., Shih, C. W.: Computer Simulation Of Impression Creep Using The 
Hyperbolic Sine Stress Law, Materials Science And Engineering, V201, Oct 1995 

22. Lee, S., Jeng, S. M., Yang, J. -M., Modeling and Simulation of the Effect of Fiber Breakage 
On Creep Behavior Of Fiber-Reinforced Metal Matrix Composites, Mechanics Of Materials, 
V21, Nov 1995 

23. Loghman, A., Wahab, M. A., Creep Damage Simulation of Thick-Walled Tubes Using The 
[Theta] Projection Concept, International Journal of Pressure Vessels and Piping, V67, Jun 
1996 


NASA/CR— 2003-212351 


98 



24. Ilsup, Chung, Sun, C.T., Chang, I.Y.: Modeling Creep In Thennoplastic Composites, Journal 
Of Composite Materials, V27, 1993 

25. Jain, M.K., Iyenger, S.R.K., Jain, R.K.: Numerical Methods For Scientific and Engineering 
Computation, New Age International Limited, Third Edition, 1996 

26. Penny, R.K., Mariott, D.L.: Design For Creep , McGraw Hill, 1971 

27. Hinton, M.J., and Soden, P.D., Predicting Failure In Composite Laminates: The Background 
To The Exercise. Composite Science and Technology, 58, 1001. 1998 

28. Gotsis, P.K., Chamis, C.C., and Minnetyan, L.: Prediction Of Composite Laminate Fracture: 
Micromechanics and Progressive Fracture. Composites Science and Technology, 58, 1137. 
1998 

29. Sun, C. T. and Tao, J. X., Prediction Of Failure Envelopes and Stress/Strain Behavior Of 
Composite Laminates. Compos. Sci. Technol., 1998, V58 

30. Rotem, A., Prediction Of Laminate Failure With Rotem Failure Criteria. Compos. Sci. 
Technol., 1998, V58. 

31. Wolfe, W. E. and Butalia, T.S., A Strain-Energy Based Failure Criteria For Nonlinear 
Analysis Of Composite Laminates Subjected To Biaxial Loading. Compos. Sci. Technol., 
1998, V58. 

32. Agarwal, B.D., and Broutman, L.J.: Analysis and Performance Of Fiber Composites, 2 nd 
Edition, John Wiley and Sons Inc. 1990. 

33. Chamis, C.C.: Failure Criteria For Filamentary Composites. NASA Tn-D-5367, 1969 

34. Tsai, S.W.; Strength Theories Of Filamentary Structures In Schwartz’s Fundamental Aspects 
Of Fiber Reinforced Plastic Composites, Interscience, Ny. 1968. 

35. Hill, R.: The Mathematical Theory Of Elasticity , Oxford University Press, Oxford. 1950. 

36. Chamis, C.C., Lark, R.F., and Sinclair, J.H.: Integrated Theory For Predicting The 
Hygrothermal Response Of Advanced Composite Structural Components. ASTM STP 658. 
1978 

37. Chamis, C.C., and Hopkins, D.A.: Thermoviscoplastic Nonlinear Constitutive Relationships 
For Structural Analysis Of High Temperature Metal Matrix Composites. NASA Tm 87291. 
1985. 

38. Gibson, R. F., Principles Of Composite Material Mechanics , McGraw-Hill, 1994, Pp. 139. 

39. Chamis, C.C., and Sinclair, J.H.: Durability/Life Of Fiber Composites In 

Hygrothermalmechanical Environments. ASTM Stp-787, ASTM, Philadelphia. 1982. 

40. Minnetyan, L., Murthy, P.L.N., and Chamis, C.C.: Composite Structure Global Fracture 
Toughness Via Computational Simulation, Computers and Structures, Vol. 37, No. 2. 1990. 


NASA/CR— 2003-212351 


99 



REPORT DOCUMENTATION PAGE 

Form Approved 
OMB No. 0704-0188 

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 

1 . AGENCY USE ONLY ( Leave blank) 

2. REPORT DATE 

June 2003 

3. REPORT TYPE AND DATES COVERED 

Final Contractor Report 

4. TITLE AND SUBTITLE 

5. FUNDING NUMBERS 


Quantification of Energy Release in Composite Structures 


6. AUTHOR(S) 

Levon Minnetyan 


WBS-22-708^8-1 1 
NAG3-2393 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Clarkson University 
1 Clarkson Avenue 
Potsdam, New York 13676 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-13939 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA CR— 2003-212351 


11. SUPPLEMENTARY NOTES 


Project Manager, Christos C. Chamis, Research and Technology Directorate, NASA Glenn Research Center, 
organization code 5000, 216-433-3252. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified - Unlimited 
Subject Category: 39 

Available electronically at http://gltrs.grc.nasa.gov 

This publication is available from the NASA Center for AeroSpace Information, 301-621-0390. 


Distribution: Nonstandard 


12b. DISTRIBUTION CODE 


13. ABSTRACT (Maximum 200 words) 

Energy release rate is usually suggested as a quantifier for assessing structural damage tolerance. Computational 
prediction of energy release rate is based on composite mechanics with micro-stress level damage assessment, finite 
element structural analysis and damage progression tracking modules. This report examines several issues associated 
with energy release rates in composite structures as follows: Chapter I demonstrates computational simulation of an 
adhesively bonded composite joint and validates the computed energy release rates by comparison with acoustic 
emission signals in the overall sense. Chapter II investigates the effect of crack plane orientation with respect to fiber 
direction on the energy release rates. Chapter III quantifies the effects of contiguous constraint plies on the residual 
stiffness of a 90° ply subjected to transverse tensile fractures. Chapter IV compares ICAN and ICAN/JAVA solutions 
of composites. Chapter V examines the effects of composite structural geometry and boundary conditions on damage 
progression characteristics. 


14. SUBJECT TERMS 

Damage tolerance; Computational methods; Computer codes 


15. NUMBER OF PAGES 

110 


16. PRICE CODE 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 


